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ABSTRACT 

We present the average metallicity and star-formation rate of Lya emitters (LAEs) measured from 
our large-area survey with three narrow-band (NB) fihers covering the Lya, [O ii] A3727, and IIa-|-[N ii] 
hnes of LAEs at z — 2.2. We select 919 z = 2.2 LAEs from Subaru/Suprime-Cam NB data in 
conjunction with Magellan/IMACS spectroscopy. Of these LAEs, 561 and 105 are observed with 
KPNO/NEWFIRM near-infrared NB filters whose central wavelengths are matched to redshifted 
[On] and Ha nebular lines, respectively. By stacking the near-infrared images of the LAEs, we 
successfully obtain average nebular-line fluxes of LAEs, the majority of which are too faint to be 
identified individually by narrow-band imaging or deep spectroscopy. The stacked object has an Ha 
luminosity of 1.7 x 10'*^ergs~^ corresponding to a star formation rate (SFR) of 14M0yr~^. We 
place, for the first time, a firm lower limit to the average metallicity of LAEs oi Z > 0.09 Zq {2a) 
based on the [O ii]/(Ha+[N ii]) index together with photo- ionization models and empirical relations. 
This lower limit of metallicity rules out the hypothesis that LAEs, so far observed at z ~ 2, are 
extremely metal poor (Z < 2 x 10^^ Zq) galaxies at the 4a level. This limit is higher than a simple 
extrapolation of the observed mass-metallicity relation of z ~ 2 UV-selected galaxies toward lower 
masses (5 x 10^ Mq), but roughly consistent with a recently proposed fundamental mass-metallicity 
relation when the LAEs' relatively low SFR is taken into account. The Ha and Lya luminosities of 
our NB-selected LAEs indicate that the escape fraction of Lya photons is ~ 12 — 30%, much higher 
than the values derived for other galaxy populations at z '-^ 2. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: formation — galaxies: 
high-redshift — galaxies: star formation 



1. INTRODUCTION 

Galaxy mass is thought to be a fundamental quan- 
tity which governs the evolution of galaxies. De- 
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tailed observations of present-day galaxies show that 
various properties of galaxies such as star formation 
rate (SFR ) and gas phase metallicity correlate with 
mass (e.g., !Brinchm ann et al.ll20Q4l:lTremonti et"aLll2004l: 
I Ellison et a l. 2008). Theoretically, galaxy formation 
models based on ACDM cosmology predict that galaxies 
grow through subsequent mergings of lower-mass objects 
and that galaxy properties are largely determined by 
their masses through mass-dependent processes at work 
in the evolution o f galaxies (e.g.,[ Blumcnthal ct al. 1985 
iDavis et all 119851 : iBardeen et al.ril986i) . Therefore, ob- 
servations of the mass-dependence of galaxy properties 
back in cosmic time are crucial to understanding how 
galaxies evolve to acquire present-day properties. 

Past observations have revealed that some physi- 
cal quantities like SFR and metallicity correlate with 
mass also at high-z, although the correlation seems 
to ev olve with redshift relative to the local galaxies 
(e.g.. [Reddv et al.l [20061: iHavashi et"all 120091: lErb et al.l 
l2006ailMaiolino et al.l l2008: Ma nnucci et al.ll2009[ ). Thev 
demonstrated the importance of mass-dependent effects 
on galaxy evolution. However, high-z samples are bi- 
ased for high-mass galaxies (> 10^ -^0) because most of 
them are based on continuum selected sampl es, such as 
Lyma n-break galaxies (LBGs ; e.g., [S teidcl fc HamiltonI 
Il992f) and BzK galaxies (e.g.. iDaddi et al.,.200ir 

Lya emitters (LAEs), galaxies commonly observed at 
high redshifts with strong Lya emission, are likely to 
be low-mass, young galaxies as suggested from their 
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TABLE 1 

Summary of NB387 Imaging Data in the SXDS field 



Field Name 


Exp. Time 


PSF 


Area 


mi 


im 


Date 




(1) 






(2) 


(3) 


(4) 




SXDS-C 


3.20 


9 




0.88 


587 [41] 


25.7 


25.3 


2009 Dec 14 - 16 


SXDS-N 


2.50 


5 




0.70 


409 [159] 


25.6 


25.2 


2009 Dec 16 


SXDS-S 


2.50 


5 




0.85 


775 [344] 


25.7 


25.3 


2009 Dec 16 


SXDS-E(5) 


3.33 [10] 


1.95 








2009 Dec 19, 20 


SXDS-W 


1.83 [5] 


1.23 


232 [122] 


25.1 


24.7 


2009 Dec 16, 19 



^"'"■^ Total exposure times (hour). The value in square brackets shows the numbers of exposures that are combined. 
(2) FVVHM of PSFs of the stacked image that are registered with broadband images (arcscc). 

EfTeetive area that is used for the selection (arcmin'^). The value in square brackets shows the area with \ow-S/N . 
^'^^ Limiting magnitude derived from 5(T sky noise in a 2" diameter aperture. Note that murn in the SXDS-W is defined 
in a 2. 5 diameter aperture. The magnitude in square brackets shows the limiting magnitudes in the \ow-S/N regions. 
(5) The SXDS-E has large PSF size and is not included in the following analysis. So, we do not report the and 
covered area in this field. 



small sizes, faint continua, and low stellar masses in- 
ferred from spectral energy distribution (SED) fitting 
(< 10^M^: IGawiser et all [2001^ [2Q07t IFinkeM^^ al l 
200l [20081 [20091: iNilsson et al.ll20Q7l: iPirzkal et al."20Q7^ 



Lai et al. lOno et aL I I2010al lbl: Yuma ct al. 2010) 



Since they can be efficiently detected by narrow-band 
imaging, LAEs are a useful probe to investigate low-mass 
galaxies in the early stages of galaxy evolution. Further- 
more, low-mass galaxies at high redshifts such as LAEs 
are especially interesting since they are likely to be build- 
ing blocks of massive galaxies seen in later epochs. 

Thanks to the remarkable progress in observations of 
LAEs, our kno wledge of their prop e rties is rapidly accu- 
mula ti ng f^e.g.. | Cowie fc HulflOO a 'Malhotra & Rhoadsl 
2002t lOuchi et al.l 1 20031: llVIal hotra & Rhoads 2004; 
Gawiser et al.l 120061: 'Kashikawa 2006; Shiniasaku ct alj 



200l" 



Gronwall et al. 2007; Ouchi ct al. 2008|_ 
Nilsson ct al. 2009; Blanc et al.1 12011: Guait a et al 



2010: Haves et al..,2010: .On 



201(1: 



o et al.l 

mm 



.Finkelstein et al.l 

Nilsson et al.l 1201 ID . However, in almost all the ob 



2010a.b.: .Ouchi et al 
Guaita et al.l 120111 : 



servations, physical properties of LAEs including SFR 
and metallicity have been generally estimated by 
the SED fitting of broadband photometry. This is 
in contrast to massive galaxies at similar redshifts. 
While a large number of massive, continuum-selected 
galaxies have now direct measuren ients of SFR and 
metallicity from nebular lines (e.g., 'Reddv et al.l 120061: 
Erb et al. 2006b; Maiolino et al. 2008; Ma nnuc ci et al l 
20091: iHavashi et al1l2009inYoshikawa et al.ll2010l) . there 



are only a few LAEs with such direct meas urements 
(jMcLinden et al.ll2011t iFinkelstein et al.ll2011| ). One of 

the major reasons for the paucity of the direct mea- 
surements is that well-studied LAEs are located at very 
high redshift (3 < z < 7), where (rest- frame optical) 
nebular lines redshift into infrared wavelengths that 
cannot be observed from the ground. However, nebular 
lines of bright LAEs at moderate redshifts (z = 2 — 3) 
have i iow been measured thr o ugh recent LAE s urvey s 
(e.g., iMcLinden et al.l 120111: IFinkelstein et al.l 1201 ll ). 
It is known that SFRs derived from SED fitting are 
dependent on the star formation history assumed, and 
they can vary by an order of ma gnitude among different 
histories (e.g., lOno et ani2010al ). It is also known that 
SED fitting cannot strongly constra in metallicities du e 
to the degeneracy with age (e.g., lOno et all l2010aD . 
Observations of nebular lines for a large number of 



LAEs are essential to extending star formation rate 
and metallicity measurements toward low-mass galaxies 
below ~ 10^ M©, so that the mass-dependencies of SFR 
and metallicity can be compared with their present-day 
homologue over a full mass range. 

We are conducting an imaging survey of z ~ 2.2 
LAEs in several fields on the sky using three narrow- 
band filters described below. This redshift is unique be- 
cause [O ii]A3727 and Ha lines fall into wavelength ranges 
where OH-airglow is very weak, thus enabling one to 
study SFRs and metallicities of LAEs using these lines 
from the ground. We developed a new narrow-band filter, 
NB387, with a central wavelength and FWHM of 3870 A 
and 94 A, respectively, to select LAEs over z — 2.14 - 
2.22. [On] and Ha lines in this redshift range are then 
observed through near-infrared (NIR) narrow-band fil- 
ters, NB118 (Ac = 11866 A, FWHM= 111 A) and NB209 
(Ac = 20958 A, FWHM= 205 A), respectively, devel- 
oped by the NewHa Survey (Lee et al. in preparation). 
In this paper, we present the results from data of the 
Suha.r\i/ XMM-Newton Deep Survey field, which are the 
first results of this triple-narrowband survey. 

While our NB118 and NB209 imaging in the SXDS is 
not deep enough to detect [O ii] and Ha lines for individ- 
ual objects except for very luminous ones, we successfully 
detect these emission lines in the stacked images of more 
than 100 LAEs. Wc then measure their fluxes to derive 
the SFR and metallicity, and place the constraints on the 
average SFR and metallicity of a typical LAE at z = 2.2 
for the first time. We use these measurements to discuss 
the SFR and metallicity of low-mass (< 10^ Mq) galaxies 
at z ~ 2. 

Since Lya photons produced in a galaxy are expected 
to be easily absorbed by dust in the interstellar medium 
(ISM) of the galaxy during repeated resonant scatterings 
by neutral hydrogen gas, a puzzle is why LAEs have such 
strong Lya emission. There are three possible answers to 
the p uzzle. The first is th at LAEs are primordial galaxies 
(e.g..lSca nnapieco et al.l [2003). If population HI stars are 
formed following a top-heavy initial mass function, hard 
UV photons would be radiated in a short time scale, and 
strong Lya lines would be observed. In this case, LAEs 
are young, and extremely metal poor galaxies. This will 
be tested by the metallicity we estimate. The second 
is th at the ISM has a clumpy geometry (e.g., iNeufeldl 
Il991h . where the Lya photons are scattered at the sur- 
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Fig. 1.— Relative response curves of the triple narrow-band filters, NB387 at Subaru and NBf 18 and NB209 at KPNO/NEWFIRM (red), 
superimposed on those of CFHT «*, Subaru B, V, ij, i', z', and UKIRT/WFCAM J, ft' (black). The responses include the throughputs 
of the instrument and the telescope as well as atmospheric absorption. The inner panels are a zoom in around NB387, NB118, and NB209, 
whose upper x-axis shows the redshift of Lyo, [On], and Ho, respectively. In the NB118 and NB209 zoom-in panels, the dashed curves 
indicate the response curves of NB118 and NB209 at the corner of the filter (see text for details; Equation l[TJ), while the dot-dashed lines 
are re-scaled response curves of NB387 along the wavelength axis to sample [O ll] and Ho after correction for a velocity offset between Lyo 
and nebular lines of 400kms~^. 

ered by deep Suprime-Cam broadband data in five point- 
ings with small overlaps; these five 'sub-fields' are named 
SXDS-C, -N, -S, -E, and -W , respectively, after the ir rel- 
ative positions on the sky ()Furusawa et al.l 1200 81) . We 
acquired NB387 imaging for all the five sub-fields (Ta- 
ble [T]). We do not, however, use the data of SXDS-E in 
this study because of bad seeing. For photometric cali- 
bration, we observed spectropho tometric s tandard stars 
Feige34, LDS749B, and G93-48 (|Okelll990 '). Each stan- 
dard star was observed more than twice under photomet- 
ric condition with airmasses of 1.1 — 1.3. 

We used the Su prime-Cam Deep Field Reduction 
package fSDFRED: lYaei et al.l [2Q0a lOuchi et al.ll2004H 



faces of the clumps and thus are not heavily absorbed 
before escaping from the galaxy. T he third is the out- 
flow of the ISM (e.g.. lKunth et al.1 11998) . Lya photons 
that are scattered at the far side of expanding ISM can be 
Doppler shifted to have redder wavelengths, and escape 
from galaxies without being heavily absorbed by neutral 
hydrogen gas. In order to investigate these possibilities, 
we estimate the Lya escape fraction, /^^" , for our LAEs 
by comparing the observed Lya luminosity with the in- 
trinsic Lya luminosity predicted from the dust-corrected 
Ha luminosity. 

This paper is organized as follows. We describe the 
data in fj2j The sample of z = 2.2 LAEs is constructed 
in fJH where the results of optical spectroscopy are also 
shown. In 21 detect [On] and Ha emission in the 
stacked LAEs and calculate their equivalent widths. Ob- 
jects with individual detections of these lines are also 
briefly mentioned. Results of SED fitting of the stacked 
LAEs are briefly described in S|5] In S|6l we derive SFR, 
metallicity, and /^^" for our LAEs, and discuss their 
implications. Conclusions are given in SJT] Throughout 
this p aper, magnitudes are given in the AB system (Ok^ 
119741 ) , and we assume a standard ACDM cosmology with 
(rJ™,rjA,i?o) = (0.3,0.7, 70kms-iMpc-^) 



2. IMAGING DATA 

In this section we describe the optical and NIR data in 
the SXDS field used in our analysis. 

2.1. NB387 Images 

We carried out NB387 imagin g observations o f the 
SXDS with Subaru/Suprime-Cam (iMivazaki eraIll2002D 
on 2009 December 14-16 and 19-20. Tabic [1] summarizes 
the details of the observations. The SXDS field is cov- 



to reduce the NB387 data. The data reduction pro- 
cess included bias subtraction, flat fielding, distortion 
correction, cosmic ray rejection, sky subtraction, bad 
pixel/satellite trail masking, image shifting, and stack- 
ing. For cosmic r ay rejection, we used LA. COSMIC 
(|van Dokkum|[200ll ). After the stacking process, our im- 
ages were registered with the archival broadband images 
( §2.2|) using bright stellar objects commonly detected in 
the NB387 and the broadband images. 

The PSF sizes of the registered images for the four 
sub-fields are 0'.'70 — 1'.'2. The 5a detection limits in a 
2" diameter aperture are 25.1 — 25.7 mag except on the 
edges of the images where signal-to-noise ratios (S/N) 
were significantly lower due to dithering. However, we 
include those low S/N regions to increase the number 
of LAEs, setting brighter limiting magnitudes according 
to the S/N ratios. The limiti ng magnitudes are est i- 
mated in the same manner as in iFurusawa et al.l ()2008l) ; 
we spread 5, 000 2'.'0 diameter apertures over the en- 
tire image randomly after masking detected objects, and 
measure their photon counts. We then fit the negative 
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TABLE 2 

Summary of Optical Broadband Imaging Data in the SXDS field 



Band Observatory Field Name PSF miim Reference 

(1) (2) (3) 



u*W CFHT SXDS-C,N,S,E,W 0.85 26.9 (a) 

B Subaru SXDS-C,N,S,E,W 0.78 - 0.84 27.5 - 27.8 (b) 

V Subaru SXDS-C,N,S,E,W 0.72 -0.82 27.1 -27.2 (b) 

R Subaru SXDS-C,N,S,E,W 0.74-0.82 27.0-27.2 (b) 

i' Subaru SXDS-C,N,S,E,W 0.68-0.82 26.9-27.1 (b) 

z' Subaru SXDS-C,N,S,E,W 0.70-0.76 25.8-26.1 (b) 

(1) pgp si2c is defined as a FWHM of point sources (arcscc). 

(2) The limiting magnitude (Sct) estimated by 2" diameter random aperture photometry, 
(a) Foucaud et al. in preparation, fb) |Furusawa et al.| i200g|) 

^ ^ The u* image covers the same area of the UKIDSS/UDS project i|L awrence et al.ll2007l) which corresponds to 
about 77% at the SXDS-C,-N,-S,-E, and -W. 



Relative R.A. [Mpc(comoving)] 



-20 -10 10 20 




R.A. [deg] 

Fig. 2. — Areas covered by respective imaging data in the SXDS. 
Lines in different colors outline areas covered in different pass- 
bands: blue: NB387; magenta: u* and J, H, K; orange: NB118; 
red: NB209; black: optical {_B, _R, i', 2'). The orange-shaded and 
red-shaded regions correspond to the NB118 sub-region and the 
NB209 sub-region, respectively. Note that the NB209 sub-region 
is embedded in the the NB118 sub- region (see Js}. 

part of the histogram of the counts with Gaussian, and 
regard its a as the la sky fluctuations of the image for 
2'.'0 diameter apertures. The limiting magnitudes which 
appear in the following subsections are also estimated in 
the same manner. 

We infer the errors in photometric zero points of our 
NB387 images from a comparison of colors (u*— NB387 
or B— NB387) of stellar objects in the images wit h those 
of 175 Galactic stars calculated from Gunn & S trvkerl 
(|1983[ )'s spectrophotometric atlas. The estimated errors 
are less than 0.05 mag, which are small enough for our 
study. 

2.2. Optical Broadband Images 

The optical broadband data are required not only 
for selecting LAEs but also for performing SED fitting 
of the stacked object (Y. Ono et al. in preparation; 
briefly mentioned in Sj5]). We use the publicly available 



B, V , R, i', and z' data take n with Subaru/Sup rime- 
Cam by the SXDS project (iF urusawa et all [2008( 1 . and 
the u* data taken with GFHT/MegaCam as part of 
the UKIDSS/UDS project (Foucaud et al. in prepa- 
ration). Table [2] summarizes the details of the optical 
broadband data. The passbands of u* and B, which 
are used to select LAEs as off-bands of NB387, are 
shown in Figure [H We register the u* image with the 
Suprime-Cam B image using common bright stars. The 
5(7 limiting magnitudes on a 2" diameter aperture are: 
26.9 (u*), 27.5-27.8(5), 27.1-27.2 (V"), 27.0-27.2 (i?), 
26.9-27.1 (i'), 25.8-26.1 (z') mag. The u* image covers 
77 % of the Suprime-Cam field, nearly t he same area as 
that of the UKIDSS/UDS JHK images (|Lawrence et all 
[20f)7l see §2.4p Figure [2] illustrates the sky coverages of 
our NB387 imaging, the Suprime-Cam broadband imag- 
ing, and the u* imaging. 

For each sub-field, PSF sizes were matched to the worst 
one among the u*BVRi'z' and NB387 images with the 
IRAF task GAUSS. The resulting PSF FWHMs of the 
images in the SXDS-C,N,S,W are O'/SS (SXDS-C), 0'.'85 
(N), 0'.'85 (S), and T.'23 (W). 

2.3. NB118 and NB209 Images 

The SXDS field has been partly imaged in 
the NIR NB118 and NB209 narrowbands with 
KPNO/NEWFIRM by the NewHa Survey (Lee et al. 
in preparatio n). More details on the NB118 observations 
are given in 'L v et al] ()201lD . The regions imaged are 
shown in Figure [21 The response curves of these two 
narrowband filters are shown in Figure [1] with those of 
the WFCAM J, H, and K bands (j fej) . Exposure times 
were 8.47 - 12.67 hr for NB118 and 11.75 hr for NB209. 
Both images were registered to the Suprime-Cam z' band 
images using common bright starj^. The ha limiting 
magnitudes in a 2" diameter aperture are estimated to be 
~ 23.6 mag for the NB118 and ~ 22.6 mag for the NB209 
[^ . Errors in photometric zero points of the NB118 and 
NB209 images are less than 0.05 mag, which are inferred 
following the same manner as the NB387 images. 

In NEWFIRM the incident angle to the filter surface 
is not exactly normal but varies as a function of the dis- 
tance from the field center (Lee et al. in preparation; 

in order to run SExtractor in double-image mode (see ^3.11 1. 
^"^ Our effective limiting magnitudes of NB118 and NB209 are 
0.2 — 0.3 mag deeper than the magnitudes listed here owing to 
smoothing procedures (see i|4.1ll. T he limiting magnitudes after 
the smoothing are described in ii4.2l 
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u*-NB387 

Fig. 3.— _B-NB387 vs. m'-NBSS? plane for selection of LAEs 
at 2 ~ 2.2. The solid lines in various colors indicate tracks of model 
galaxies redshifted from to 3.5 with and without Lya emission: 
purple: simple stellar population with an age of 0.03 Gyr; green, 
cyan, and blue: simple stellar population with Lyo emission of 
EWrost = 200,100, and 40 A, respectively. The symbols on the 
tracks of model galaxies with hya emission correspond to z = 
2.14 (filled triangles), 2.16 (filled squares), 2.18 (filled circles), 2.20 
(open squares), and 2.22 (open triangles). The red and orange 
solid lines represent tracks of elliptical and spiral galaxies from the 
SWIRE template library (PoUetta et al. 2007), respectively. The 
yellow stars show Galactic stars from Gun n fc Strvker, (,1983). The 
tilted arrow indi cates the redd ening effect in the case of E{B — 
V) = 1.0 iCalzetti et al.ll200a) . The objects that are located in 
the area enclosed by the solid black lines are considered to be LAE 
candidates. 



see also lTanaka eFaLl (|20T1 for Subaru/MOIRCS). Ac- 
cordingly, the central wavelength of NB118 and NB209 
also varies over the FoV, since both filters are interfer- 
ence filters. The central wavelength at an incident angle 
9 (angle from normal incidence) is given by 



X{0) = Ao cos(6'/n) 



(1) 



where Aq is the wavelength at normal incidence and n 
is the index of refraction of the material on which inter- 
ference film is coated. We adopt n = 1.50 for NB118 
and n — 1.49 for NB209 (Lee et al. in preparation). 
At the corner of the filter 9 has the maximum value 
of 13.5 degree. Thus, the maximum wavelength shifts 
are ~ -150 A and ~ -260 A for NB118 and NB209, re- 
spectively. The passband of NB387 is broad enough to 
cover Lya lines at the redshifts of [O ii] and Ha lines 
corresponding to these shifted wavelengths (see also in- 
ner panels of Figure [ij. However, as seen in M.1.21 we 
have to take into account the passband shift over the FoV 
when calculating the equivalent widths of [O ii] and Ha 
lines of stacked objects. 

2.4. NIR Broadband Images 

The UKIDSS/UDS project provides deep J, H, an d K 
images of UKIRT/WFCAM ([Lawrence et all [2001 . In 
this paper, we use the data release 8 (DR8) images cur- 
rently available to the UKIDSS consortium. The J and 



K images are used as off-bands of NB118 and NB209 to 
detect [On] and Ha(-|-[Nii]) emission, respectively, while 
the H image is used for SED fitting of the stacked ob- 
jects. We register these three images with the Suprime- 
Cam z' band images in the same manner as for NB118 
and NB209. The 5cr limiting magnitudes over a 2" di- 
ameter aperture are estimated to be 24.8, 24.1, and 24.6 
in the J, H, and K bands, respectively. J, H, and K 
response curves are shown in Figure [1] 

3. PHOTOMETRIC SAMPLES OF LAES AT Z = 2.2 

We select LAEs in a 2, 003 arcmin^ region which is 
covered by all the three passbands for selecting LAEs: 
NB387, u* , and B (see Figure[2|). In this paper, however, 
we use only LAEs in a sub-region of 1, 283 arcmin^ which 
is also covered by NB118. In this 'NB118 sub-region', a 
353 arcmin^ region is covered by NB209 as well ('NB209 
sub- region'). Note that the J, H, K data are available for 
the NB118 sub-region. We use LAEs in the NB118 sub- 
region to derive the typical [O ii] luminosity of LAEs, and 
those in the NB209 sub-region to derive the typical Ha 
luminosity of LAEs. The typical [O ii] luminosity is also 
derived for LAEs in the NB209 sub-region to estimate 
the metallicity in combination with the Ha luminosity. 
In this section, we describe the construction of the LAE 
sample for the entire 2,003 arcmin^ region. 

3.1. Object Detection and Candidate Selection 

We use the SExtractor software ([Bertin fc Arnout^ 
Il996( ) to perform source detection and photometry. We 
run SExtractor in double-image mode, with the NB387 
image used as the detection image. We identify sources 
with 5 adjoining pixels and brightness above > 2a of the 
sky background. The diameter to define aperture mag- 
nitudes is set to 2'.'5 for the SXDS-W and 2'.'0 for the 
other sub-fields. We use aperture magnitudes to calcu- 
late colors, and adopt MAG_AUTO for the total magnitude. 
All magnitudes are corrected for Galact ic extinction of 
E{B -V) = 0.020 (ISchlegel et all [19981) . The NB387- 
detection catalog contains 42, 995 objects with aperture 
magnitudes brighter than the 5a sky noise. 

We select LAE candidates on the u*-NB387 vs. 
i?— NB387 color plane (Figure [3]). In this figure, colors 
of model galaxies and Galactic stars are plotted in order 
to define the selection criteria for LAEs. The tracks indi- 
cate the colors of model galaxies redshifted from 0.00 to 
3.50 with a step of Az = 0.01. We assume the inter galac - 
tic medium (IGM) attenuation model of Ma daul ([19951 ) . 
Base on Figures [3| and [5l we define the color criteria of 
z - 2.2 LAEs as: 



NB387 > 0.5 kkB- NB387 > 0.2, 



(2) 



which select LAEs with EWrest > 30 A. The 2a pho- 
tometric errors in m*— NB387 for the faintest objects 
(NB387=25.7) in our NB387-detected catalog are ~ 0.5 
mag. Thus, the criterion of u*— NB387> 0.5 ensures that 
the contamination fraction in our LAE sample due to 
photometric errors is sufficiently low. We use 2a limit- 
ing magnitudes instead when an object is not detected 
in u* or B at 2a level. The selection criteria require 
NB387 magnitude significantly brighter than both u* and 
B magnitudes, which results in small number of non- 
emitters in the sample. 
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Fig. 4. — Distribution of NB387-detected objects in the 
B-NB387 vs. «*-NB387 plane. The black dots indicate all the 
detected objects, while the red filled circles show LAE candidates 
after removing spurious objects and interlopers. The blue open 
squares with errorbars show spectroscopically confirmed LAEs. For 
the purpose of display, objects whose «*— NB387 colors exceed 2.4 
are plotted at «*-NB387 = 2.4. 

Using Equation © , we identify 1 , 044 LAE candidates 
in the NB387 detected catalog. These candidates are 
contaminated by spurious objects and foreground and 
background interlopers. We remove those contaminants 
by the procedures described in the next two subsections. 

3.2. Rejection of Spurious Objects 

Given the limited number of ditherings in our NB387 
imaging (~ 4 — 5), a clipped- mean stacking fails to com- 
pletely remove the remaining cosmic rays which survived 
LA. COSMIC rejection on individual dithered images. 
Since cosmic rays have steeper light profiles than the 
PSF, we removed 21 sources with FWHMs significantly 
smaller than the PSF. We then perform visual inspection 
on all the remaining objects, and eliminate 90 obvious 
spurious sources such as ghosts due to bright stars and 
bad pixels. 

3.3. Identification of Interlopers 

In addition to LAEs, other emission line objects, such 
as [On] emitters at z ~ 0.04, MgiiA2798 emitters at 
z ~ 0.4, and CivA1550 emitters at z ~ 1.5, may be 
included in our sample. Our survey area, 2, 003 arcmin^, 
corresponds to 400 Mpc'^ for [O ii] emitters, which is two 
orders of magnitude smaller than the volume sampled 
by Lya emitters (48, 100 Mpc^). The number of [On] 
emitters is therefore expected to be small. 

To remove [O li] emitters from our sample, we use 
the Galaxy Evolution Explorer (GALEX) NUV (A^ = 
2267 A, FWHM = 616 A) and FUV (Ac = 1516 A, 
FWHM = 269 A) data. Real LAEs in our sample will 
be invisible in these data, since these two passbands are 
located shortward of the Lyman break at z ^ 2.2. Thus, 
objects visible in either of the two GALEX band data are 




20 21 22 23 24 25 



NB3B7 

Fig. 5. — Distribution of NB387-detected objects in the 
n*-NB387 vs. NB387 plane. The black dots indicate all the de- 
tected objects, while the red filled circles show LAE candidates. 
The blue open squares with errorbars show spectroscopically con- 
firmed LAEs. For the purpose of display, objects whose u*— NB387 
colors exceed 2.4 are plotted at «*— NB387= 2.4. The horizontal 
solid line shows the selection threshold of «*— NB387 and the blue 
curve indicates the 2a photometric error in u*— NB387 for sources 
with «*— NB387 = 0.05, which is the average «*— NB387 color of 
all the objects. The right y axis shows the rest-frame Lya equiva- 
lent width of z = 2.18 LAEs with «*— NB387 color corresponding 
to the left y axis. 

likely to be [O n] emitters. Although some studies have 
shown that ionizing photons are more likely to escape 
from Lya selected ga l axies than from UV -selected galax- 
ies fe.g.. Ilwatal 120091: INestor et al.ll20lil ). the estimated 
UV-to- Ly-continuum flux density ratio is > 2 even for 
LAEs ([Nestor et al.|[2QTTI ). therefore the Ly-continuum 
oi z ^ 2 LAEs is expected to be fainter than the UV 
continuum by ^ 1 mag or more. Since the detection 
Hmit of the GALEX images is ~ 24 (3cr) both in NUV 
and FUV, LAEs fainter than V ^ 23 should be invis- 
ible in the GALEX images. We find 12 objects which 
have a counterpart in either of the NUV or FUV image 
within 3" from the NB387 position. Among them, 4 are 
fainter than V = 23, thus to be inferred to be interlop- 
ers. The remaining eight are brighter than V = 23, and 
their GALEX magnitudes are also bright enough to be 
consistent with them being interlopers. Indeed, all eight 
are also detected as an X-ray or radio source as described 
in the next paragraph. We thus remove these 12 objects 
from the sample. 

For Mgii and Civ emitters, we use X-ray and radio 
imaging data, since both emitters selected by Equation 
© should have large EWs (EW^cst > 30 A) and thus 
they are likely to be AGNs. For X-ray data, we use 
the XMM-New t on 0.2 - 10.0 keV band source catalog 
by lUeda et all (|2008f ). For radio data, we use the 
Very Large A r ray (V LA) 1.4 GHz source catalog by 
iSimpson et all ()2OO60 . After removing some confused 
objects by visual inspection, we find 11 (1) LAE candi- 
dates with an X-ray (radio) counterpart. The number 
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TABLE 3 

Properties of the spectroscopic sample 



ID 


R.A.(i) 




mag(NB387) 


-^obs 


z 


flag 








(2) 


(3) 


(4) 


(5) 


NB387-C-04640 


02:18:48.968 


-05:09:00.32 


23.07 ± 0.02 


3862.83 


2.1767 


b 


NB387-C-07673 


02:18:56.532 


-05:05:48.41 


24.98 ± 0.07 


3853.55 


2.1690 


c 


NB387-C-08099 


02:19:05.729 


-05:05:23.86 


24.79 ± 0.05 


3889.02 


2.1982 


c 


NB387-C-08204 


02:18:57.385 


-05:05:18.82 


24.46 ± 0.05 


3895.63 


2.2036 


b 


NB387-C-08321 


02:19:05.279 


-05:05:11.22 


24.89 ± 0.06 


3891.62 


2.2003 


b 


NB387-C-09219 


02:19:02.396 


-05:04:19.27 


24.37 ± 0.04 


3890.03 


2.1990 


c 


NB387-C-09951 


02:18:50.038 


-05:03:34.09 


24.38 ± 0.04 


3901.78 


2.2087 


b 


NB387-C-11135 


02:18:37.381 


-05:02:24.61 


23.86 ± 0.03 


3882.13 


2.1925 


b 


NB387-C-12596 


02:18:55.071 


-05:00:58.82 


24.83 ± 0.06 


3886.29 


2.1960 


c 


NB387-C-16564 


02:19:09.542 


-04:57:13.32 


22.88 ± 0.03 


3861.44 


2.1755 


a 



Coordinates arc in J2000. 

NB387 magnitudes with 2" diameter apertures. 
^■^^ Central wavelengths of observed Lya lines (A). We perform a Gaussian fitting to each detected line using 
MPFIT, and derive the central wavelengths. 

'•'^^ Redshifts of Lya lines estimated from the central wavelengths. 



(5) 



Flags of reliability of the lines by inspecting the 2D spectra by eye: "a" secure, "b" likely, "e" tentative. 
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Fig. 6. — IMACS spectra of ten confirmed LAEs. For each object, the top panel shows the 2D spectrum, while the bottom panel showing 
the ID spectrum (blue line) and the sky background (gray shaded area), both of which have been smoothed with a 3 pixel boxcar filter 
and arbitrarily normalized. 



of likely Mgii or Civ emitters is thus 12. Note that 10 
out of the 12 are also detected in the GALEX data. 

After removing spurious objects and obvious [O ii] and 
Civ emitters, we have 919 (= 1, 044 - (21 -t- 90) - (12 + 
12 — 10)) LAE candidates. Among them, 561 are in the 
NB118 sub-region and 105 in the NB209 sub-region. We 
plot all the candidates on the u*-NB387 vs. B-NB387 
two color plane in Figure |4] and on the it*— NB387 vs. 



NB387 color-magnitude plane in Figure [S] Figure |3] and 
FigureSlshow that the candidates (red circles) are placed 
in the isolated region away from the locus of other galax- 
ies and Galactic stars. 

3.4. Follow-up Spectroscopy Data 

We carried out follow-up spectroscopy of 30 objects 
selected from the whole sample {N = 919) so that they 
are distributed in wide ranges of NB387 magnitudes and 
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u*— NB387 colors. The observations were made with 
the Inam ori Magellan Areal Camera and Spectrograph 
(IMACS: [Dressier et aI]|2006D on the Magellan telescope 
using the 300 lines mm~^ grism and the WB3600 — 5700 
filter on 2010 July 9-10 under photometric conditions. 
We used the //4 camera, which has a better sensitivity 
than the //2 camera at short wavelengths. The total 
throughput of the //4 camera with the 300 lines mm~^ 
grism at 3870A is 6.2 %, while that of the //2 camera is 
2.6 %. The on-source exposure time was 13, 700 seconds, 
with a seeing size of 0'.'48 - 0'.'63. We chose a 0'.'8 slit 
width, which gives a resolving power of _R ~ 700 around 
4000 A. The COSMOS pipeline was used for data reduc- 
tion. 

By inspecting the reduced spectra by eye, we detected 
an emission line around 3870 A for ten objects, while 
the remaining 20 had no visible emission line. The 
main reason for this low detection rate is the bright 
limiting flux of our observation due to the low sensi- 
tivity below 4000 A and the relatively short exposure 
time. Indeed, all ten objects with line emission are 
brighter than NB387 = 25.0, and the success rate lim- 
ited to NB387 < 25.0 (13 objects in total) is found to 
be 10/13 = 77%. Among the remaining three bright 
objects, two have NB387~ 25.0 and may be marginally 
undetected. The other one object has NB387 = 24.4, but 
smaller Lya flux expected from its w*— NB387 compared 
with the confirmed candidates. The spectra of the ten 
confirmed LAEs are shown in Figure [HI and the NB387 
magnitudes and Lya-based redshifts are given in Table 

m 

The ten objects are not [O ii] emitters at z ~ 0.04 be- 
cause of the lack of [O iii] A5007 line at the corresponding 
wavelength, ~ 5200 A. They are not AGNs either, with 
Mgii emission at z ~ 0.4 or Civ emission at z ~ 1.5, 
because of the absence of emission lines at 3600 -5700 A 
(e.g., Mgll: [Oil] line at ~ 5150 A, Civ: HeiiA1640 line 
at ~ 4100 A and Ciii]A1909 line at ~ 4770 A) and be- 
cause we have removed AGN candidates from the sample 
in advance (see i 33.3|) . Therefore, we conclude that all ten 
objects are LAEs at z ~ 2.2. This demonstrates that the 
contamination in our LAE sample is very low, at least 
for bright objects. 

4. [On] AND Ha EMISSION LINES 

Our NIR images are not deep enough to study [On] 
and Ha properties of LAEs based on individual detec- 
tions. Indeed, only 10 objects have detection of [Oil] 
emission and only 3 have detection of Ha-|-[Nii] emission, 
as described in M.2\ We therefore carry out a stacking 
analysis of the whole LAE sample in the NB118 sub- 
region, and discuss average emission-line properties of 
z ~ 2.2 LAEs. The individually detected objects are 
discussed below and compared with the stacked objects. 

4.1. Stacking Analysis 

Stacking is done separately for the NB118 sub-region 
and the NB209 sub-region. For the NB118 sub-region, 
we stack the NB118 and J images at the positions of 561 
LAE candidates. Similarly, for the NB209 sub-region, 
the NB118, J, NB209, and K images are stacked at the 
positions of the 105 LAE candidates. Before stacking, we 
mask regions affected by a ghost, a stellar halo, and bad 



pixels, as well as regions with relatively large noise. The 
stacked image in the NB118 sub-region are used to detect 
the [O ii] flux at the highest S/N ratio, while the stacked 
image in the NB209 sub-region are used to compare [O ii] 
and Ha fluxes in a common sample. To derive the aver- 
age Lya flux, we also stack NB387 and u* images for each 
sub-region. All stackings are done by median-stacking. 
Before the stacking, we smooth the images with gaussian 
filters so that both narrowband and broadband images 
have the same PSF sizes. This enables us to measure 
the colors of the stacked objects by aperture photometry 
(i |4TT|) . The PSF size before smoothing is 1'.'6 for NB118 
and J, 1'.'2 for NB209 and K, and 1'.'2 for NB387 and u* . 
We do not remove the individually detected objects from 
the stacking analysis in order to increase the number of 
candidates for stackings. In fact, results after removing 
the individually detected objects are consistent with the 
results in Equation (|3]) owing to the median-stacking. 

The stacked NIR images are shown in Figure [T] A sig- 
nal is clearly visible in all the passbands including NB118 
and NB209. We measure the magnitudes and colors of 
the stacked images in i)4.1.11 and convert them into line 
fiuxes using Monte Carlo simulations in i i4.1.2l 

4.1.1. Photometry 

We measure aperture magnitudes of the stacked ob- 
jects using the IRAF task phot. The aperture diame- 
ter is set to be 2'.'5 — 3'.'2 (about twice the PSF size) 
to calculate colors, while larger apertures of 4'.'8 — 6'.'1 
(depending on the passband) are adopted to obtain to- 
tal magnitudes. Errors in the magnitudes ar e estimated 
in the same manner as in lOno et al.l (|2010bl ): we create 
1, 000 median-stacked sky noise images, each of which is 
made of 561 (105) randomly-selected sky noise images in 
the NB118 (NB209) sub-regions. We then fit the nega- 
tive part of histograms with a Gaussian, whose FWHMs 
are used to estimate the limiting magnitudes. The mag- 
nitudes and errors obtained are summarized in Table ID 
The uncertainties of zero points are not included in the 
errors. 

4.1.2. Equivalent Widths and Fluxes of the Lines 

The J-NB118 and i4:-NB209 colors for the stacked 
objects are calculated to be: 

J - NB118(NB118sub) = 0.43 ± 0.05 
J-NB118(NB209sub) = 0.39±0.08, (3) 
K - NB209(NB209sub) = 0.65 ±0.11 

where the color with 'NBllSsub' and 'NB209sub' in 
round brackets is derived from the stacked object in the 
NB118 sub-region and that in the NB209 sub-region, re- 
spectively. These large, positive values with the small er- 
rors indicate significant detection of the lines. In order to 
estimate the [On] and Ha-|-[Nii] equivalent widths and 
fiuxes with the best accuracy, we used Monte Carlo sim- 
ulations. If both a narrowband and a broadband (which 
brackets the narrowband) filters have ideal top-hat re- 
sponse functions, an EW of an emission line falling in 
the narrowband can be calculated as: 

_ {fx - fx) AA^AA^ 1 
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TABLE 4 

Properties of the stacked LAEs in the two sub samples 



sample 
■■■IN]W 


Line 


apcr 

(2) 


^apcr 

(3) 


N 
™-total 

(4) 


EWrest 

(5) 


Flux 
(6) 




Luminosity 
(7) 


NB118 sub 


[On] 


25.18 ± 0.02 


24.74 ±0.05 


24.57 




9.65to |g X 10- 


18 


3.54+O-i^ X 10« 


• ■ ■ [561] 


Lya 


25.83 ± 0.01 


25.05 ±0.01 


24.87 




5.76lO-;j^ X 10- 


17 


2.11+°Pj X 10« 


NB209 sub 


[On] 


24.72 ± 0.03 


24.33 ±0.07 


24.11 


96lg 


1.44lHo X 


17 


5.26+°oil X lO^i 


• ■ • [105] 


Ha+[Nii] 
Ha(8) 

Lya 


24.72 ± 0.03 
25.75 ± 0.02 


24.07 ±0.10 
25.09 ± 0.02 


23.62 
24.94 


r,7-,+242 

^'-'--104 

63tl 


2.181°^^ X 10- 
2.06l°:3i X 10- 
4.90lo;ig X 10- 


17 
17 
17 


7.98+}-^^ X 10*1 
7.55+i-l^ X 10« 
1.80+n? X 10*2 



Numbers of stacked objeets in square brackets. 
^'^^ Broadband aperture magnitudes and their Icr errors. Broadbands are J for [On], K for Hq + [Nii], and u* 



for Lya. The diameter of 
and NB387 for Lya. 



the aperture is 2. 5 — 3. 2. Aperture sizes are chosen to have twice the PSF sizes. 

^■^^ Narrowband aperture magnitudes and their la errors. Narrowbands are NB118 for [On], NB209 for Ha+[Nll] 
The diameter of the aperture is 2. 5 — 3. '2. Aperture sizes are chosen to have twice the PSF sizes. 

Narrowband total magnitudes. The diameters of the aperture are 4'/8 and 6''l for [On] and Lya line in the NB118 sub-region, 
respectively, and 5. 7, 5. 7, and 6. 1 for [On], Ha+[Nn], and Lya line in the NB209 sub-region, respectively. Aperture sizes are chosen to 
include close to 100% of the flux. 

Rest-frame equivalent width of the lines (A). For [On] and Ha-t-[Nll], EWs are estimated by Monte Carlo simulations (see i|4.1.2l in 
details), and for Lya, the colors of u*— NB387 are used for the estimates of EWs. 

Line fluxes in unit of ergs— ^ cm— ^. 
^'''^ Line luminosities in unit of ergs- ^. 

(8) The contribution from [Nn]AA6584, 6548 lines is subtracted from the observed Ha-f [Nn] luminosity using the metallicity estimated in 
" ■ {see also ilG.l.lt . 

([!]) ; see ^'2.'S\i . These issues must be considered to obtain 
correct EW values, since the objects used for stacliing 
have different redshifts (corresponding to the different 
locations within the band width of NB387) and they are 
distributed across the FoV. In general, Equation is 
correct only when a line is located where the response 
function of the narrowband peaks. In reality, however, 
a large fraction of the LAEs are expected to have [O ii] 
and IIa+[Nii] lines off the peak responses of NB118 and 
NB209 based on the re-scaled shapes of NB387 (see Fig- 
ure [Ij. This means that the EWs of stacked objects 
calculated by Equation (j4]) are always underestimated. 

We carry out IMonte Carlo simulations taking into ac- 
count of the above two points more accurately to esti- 
mate the EWs from the observed colors of the stacked 
objects. As an example, we describe below the simula- 
tions for the [On] line of the objects in the NB118 sub- 
region. The simulations for [On] and Ha+[Nii] lines in 
the NB209 sub-region are essentially the same. 

The simulations are carried out under the assumption 
that all 561 objects have identical spectra, i.e., the same 
EW and the same underlying continuum spectra. For the 
continuum spectra, we use the best-fit spectrum from the 
SED fitting to the stacked LAE (see SgO- We then vary 
EWrest over —100 A and 300 A with 1 A interval, and for 
each value we carry out a IVIonte Carlo simulation de- 
scribed by steps 1-4 below to derive the probability 
distribution of J-NB118 color for that EW. We thus 
simulate the relation between EW and J— NB118. The 
EW of the real, stacked object is calculated by substitut- 
ing the observed J— NB118 color into the relation. 
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Fig. 7. — Snapshots of the stacked LAE in J (top left) and 
NB118 (bottom left) in the NB118 sub-region, and J (top center), 
NB118 (bottom center), K (top right), and NB209 (bottom right) 
in the NB209 sub-region. Each image is 15" X 15" in size. 

where f\ is the flux density per unit wavelength, AA 
is the width of a given filter, z is redshift, and super- 
scripts N and B indicate narrowband and broadband, 
respectively. To derive this formula, we have assumed 
that the fiux density of the continuum emission is con- 
stant over the whole wavelength range. If this formula 
is used, the colors obtained above are converted into 
EWrest, NBllSsub ([Oii])= 25 A, EWrest ,NB209sub([Oll]) — 
23 A, and EWi.cst,NB209sub(Ha+[N ii])= 58 A, respec- 
tively. 

Although these values are useful as zero-order esti- 
mates, their accuracy is not sufficient for our purpose. In 
fact the assumption on the shape of the passband used 
to derive Equation ^ is over simplified in two aspects: 
the actual NB118 and NB209 passbands are not top-hat 
but rather close to a triangle shape and the central wave- 
lengths of these passbands vary over the FoV (Equation 



" We also try two independent SEDs for the simulation; 
the best-fit SED of spectroscopically confirmed BX gal axies at 
2.0 < z < 2.2 with constant star formation history (G uaita et al.l 
120111) and the best-fit spectrum from the SED fitting to the ^ = 3.1 
stacked LAE lOno et al.ll2010"al) . and find that the results are con- 
sistent with each other well within their Icr errors. 
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Fig. 8. — Relation between EWrcst and expected J— NB118 and ii'— NB209 colors from the Monte Carlo simulations: (left) EWrcst ( [O ll] ) 
vs. J-NB118 in the NB118 sub-region; (center) EWrcst ([O ll]) vs. J-NB118 in the NB209 sub-region; (right) EWrcst (Ha-f[N ll]) vs. 
_ft:-NB209 in the NB209 sub-region. The color contours indicate the probability of J-NB118 or ii'-NB209 at given EWrcst; redder colors 
mean higher possibilities as shown in the color bars in the top. The solid black lines correspond to the measured J— NB118 or A"— NB209 
colors and the dotted black lines are their Icr errors. 



1. We first generate 561 LAEs with a given EW value, 
and assign to each one of the 561 positions on the 
FoV of the real objects without duplication. Then, 
assuming that the (Lya) redshift distribution of 
our LAEs is same as the shape of the NB387 re- 
sponse function, with a peak at z = 2.18 and an 
FWHM of IS.Z — 0.075, we randomly select a red- 
shift for each of the 561 objects from this distribu- 
tion. 

2. Spectroscopic observations of high- redshift star- 
forming galaxies have found that the redshift mea- 
sured from Lya line is offset from the nebular-line 
redshift by a few hundred kms~^, due to radial ac- 
celeration of the circ umgalactic gas emitti ng the 
line ('e.g.. iPettini et a l. 2001; Stcidcl ct al. 120101 : 

iMcLinden et al.ll2011l : iFinkelstein et al..,201ll ). We 
randomly assign a redshift offset to each of the 
561 objects simulated in step 1, assuming that the 
offsets obey the dis tribution function obtained by 

ISteidel et afl (|2010[ ) for z - 2 galaxieQ Each ob- 
ject is thus given an [On] redshift. 



For each object we calculate 



and /^^^^* using 



the EW, continuum spectrum, [O ii] redshift, and 
the NB118 response function at the position of the 
object on the FoV. We then sum the 561 objects' 
// and /^^^^®, and divide the former by the latter 
to obtain the J— NB118 color of the stacked object. 

4. We repeat steps 1 - 3 500 times to obtain the prob- 
ability distribution of J— NB118 color for the given 
EW. 



Although LAEs are so far found to have systematically 
smaller velocity offsets than LBGs, ~ 15 0kms~^ for LAEs 
HMcLinden et al.ll20Tll: IFinkelstein et al.ll20Tll ') while ~ 400 km s"! 
for LB Gs USteidel et al.l 12010 ). we use the data of ISteidel et al.l 
H2010I ) since it is based on much a larger number of measurements 
(> 40 LBGs while 4 LAEs) and thus statistically more reliable. 
In any case, the wavelength shift caused by velocity offset is much 
smaller than that due to the variation of the response curve over 
the FoV (~ ISA by a velocity offset of ~ 400 km s"'^ while up to 
~ 150A by the positional variation for NB118). 



Figure [5] shows the results of the simulations, i.e., the 
distribution of J-NB118 and i4:-NB209 color as a func- 
tion of EW. In each panel, the color contour indicates the 
probability of J-NB118 or if-NB209 at given EWs, and 
the solid line and two dotted lines indicate the central 
value and the Icr lower and upper limits, respectively, of 
the observed color. We translate the observed colors and 
their errors into EWs in the following manner. First, we 
assume that for each of J-NB118 and i4:-NB209, the 
probability distribution of the true value is a Gaussian 
with its a equal to the observed Icr error in that color. 
Next, we randomly select a color following this Gaussian 
probability distribution, and assign to it an EW based 
on the EW probability distribution against color shown 
in Figure m We repeat this procedure 10,000 times. Fi- 
nally, we sort the 10, 000 EWs in ascending order, and 
regard the 50 %-tile, 16%-tile, and 84 %-tile as the central 
value, IfT lower limit, and lu upper limit, respectively. 
From the translation procedures, we obtain EWrost for 
[On] and Ha-K[Nii]: 



EW„st,NBii8sub([0II]) = imtltk 

EW„st,NB209sub([OII]) = 96t23 A. 
EWrcst,NB209sub(Ha + [Nil]) 



(5) 



071+242 X 



The errors in EW correspond to Icr. The large scatter in 
EWrest,NB209sub(IlQ!-f[N ii] ) is duc to the large photomet- 
ric error of fs:-NB209 (see Equation ^) and the large 
variance of the NB209 passband over the FoV. 

We calculate the [O ii] and Ha fluxes from the EWs ob- 
tained here combined with the total magnitudes derived 
in !j4.1.11 and convert them into luminosities. These val- 
ues are summarized in Table S) 

We estimate the Lya EW in a similar manner to what 
we explain above; we run the same Monte Carlo simula- 
tion skipping step 2, and then translate the u*— NB387 
color into the Lya EW. The Lya EWrost derived for 
the NB118 (NB209) sub-region is 86l^ (63l^) A, from 
which we obtain the observed Lya luminosity to be 



2.11 



+0.04 



X 10 



42 



-0.03 ^ ^" (l-^0_o.07 



From previous 
L(Lya) = 10*^ - 



+0.05 X ;[^q42) ergs-i (TableSD. 
narrow-band surveys, LAEs with 
ergs^^ 



10 



,44 



are referred to as typi- 
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TABLE 5 
Properties of individual objects 



ID 


R.A.(i) 


Dec.tl) 






EWroiil 


^[OTTl 


EWHa + [NII] 


-^Ha + [NII] 








(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NB387-C-06908 


02:17:57.113 


-05:06:37.06 


45+? 




977+154 
^' '-93 


d. JU_Q 22 






NB387-C-22326 


02:17:11.390 


-04:51:14.10 


54tf, 








^•^"^-0.34 


NB387-C-26967 


02:17:19.875 


-04:46:54.76 




U.MO_g 28 


22110 


9 97+0.56 
^■^'-0.65 






NB387-N-01923 


02:17:55.840 


-04:47:40.11 


1211^2 


4.87lg-i| 


127tg 


n 01+0. 18 
3-O1-0.21 






NB387-S-25047 


02:17:32.417 


-05:12:51.09 


31+3 


19.43t}:i 


18+i 








NB387-W-00165 


02:17:07.486 


-04:53:53.75 


59l« 


■J-*O_0.64 


ot; + 19 
•^''-15 


1 40+°'32 






NB387-W-00372 


02:17:04.957 


-04:45:35.62 


32i^ 




20+? 


1-631^° 






NB387-W-02225 


02:16:46.049 


-04:59:05.44 




-I 09+0.44 


10+? 


0.86l0:« 






NB387-W-04041 


02:16:28.164 


-04:45:17.72 


63ll« 




7r; + 19 
'3-16 


07+0.27 

3-»'-0.31 


40611 


6.80«;^o 


NB387-W-04492 


02:16:23.751 


-04:57:57.92 


43+11 


9 0-1+0.55 
^■•^-'--0.49 


92ti 


9 Qo+0.26 




2.08lO;« 


NB387-W-06136 


02:16:07.936 


-05:00:07.96 


40^7 


c c:n+0.94 


121^3 


9 44+0.47 
^■^^-0.50 







Coordinates are in J2000. 
^'^^ Rest-frame equivalent width of Lya in unit of A. 
^■^-^ Lya luminosity in unit of 10^^ ergs~^. 

Rcst-framc equivalent width of [O 11] in unit of A. 

[Oil] luminosity in unit of 10*^ ergs^^. 

Rcst-framc equivalent width of Ha+[Nll] in unit of A. 

Hq+[Nii] luminosity in unit of 10*^ ergs"^. 



cal LAEs (e.g., iGronwall et all [20071 : lOuchi et al.|[2008[) 

Therefore, the stacked LAEs we obtain are expected to 
have average properties of LAEs with L(Lya) > lO**^. 

4.2. Individual Objects 

Seventeen out of the 561 LAE candidates in the 
NB118 sub-region are detected in NB118 at > 5cr levels, 
among which ten have J-NB118 > 0.0 and NB118 < 
NB118(5crH indicative of the presence of [Oil] emis- 
sion. A detailed calculation using the continuum spec- 
trum from the SED fitting shows that these detection cri- 
teria correspond to EWrost([0 ii]) > 7A and /([On]) > 
1.2 X 10~^^ ergs~^ cm~^. Similarly, seven out of the 105 
objects in the NB209 sub-region are detected in NB209 
at > 5(7 levels, among which three has _fir— NB209 > 0.3 
and NB209 < NB209(5cr), indicative of the presence of 
Hq;-|-[Nii] emission at more than the 5(t level. These de- 
tection criteria correspond to EWrost(HQ;-|-[N ii]) > 22 A 
and /(Ha-h[Nii]) > 2.1 x lO-^"^ ergs-^cm-^. Two 
of the Ha-|-[Nii] detected objects (NB387-W-04041 and 
NB387-W-04492) are also detected in [On], while one 
(NB387-C-22326) is not detected in [On], probably due 
to a variation of dust extinction and the complicated ori- 
gin of the [O ii] luminosity (see ij6.1.2p . The EWs and line 
luminosities of these objects are summarized in Table [5l 

4.3. Average EWs and Luminosities of [Oil] and Ha 

Before presenting the physical properties derived from 
[O il] and Ha emission lines, we summarize the strengths 
of these lines of the stacked LAEs. 

The stacked object from the NB118 sub-region has 
EW([Oii]) = lOetJ^A and L([Oii])= 3.54t[;;}i x 
lO^^ergs"^. This EW([Oii]) is much larger than those 
obtained for typical high -z galaxies (e.g., z ^ 1 galaxies 
from the DEEP2 survev: iCooper et al.ll2006l ). 

^0 The 5cr limiting magnitude here is defined with aperture sizes 
of 3^2 diameter for NB118 and 2"5 diameter for NB209, and esti- 
mated to be ~ 22.9 and ~ 22.4, respectively. 



The stacked object from the NB209 sub-region 
has EW(Ha-l-[Nii]) = 27ll:io4A and L(Ha-|-[N ii]) = 
7.98ti]l X 10'*^ergs"\ which correspond to EW(Ha) 
= 256j:gf A and L(Ha) = 7.55l:i iJ ergs'^ after sub- 
traction of [Nil] emission (see t j6.1.ip . This EW(Ha) 
is la rger than t hose obtained for other high-z galaxies; 
e.g., lErb et aTl (2006b) performed Ha spectroscopy of 
z ^ 2 UV-selected galaxies to find the median EW(Ha) 
~ 170 A (Note that EW(Ha) of UV-selected galaxies 
from Ha spectroscopy may be biased high, since objects 
with strong er Ha e mission can be observed more easily) . 
Recently, Cowie et al. (2011) have found that the bulk 
(~ 75%) of the local LAEs have EW(Ha) > 100 A. Thus, 
it seems that a large Ha EW is a common character of 
LAEs irrespective of redshift. 

5. SED FITTING 

We perform SED fitting for the stacked objects in 
the NB118 sub-region and the NB209 sub-region to in- 
fer their stellar populations. We note that the stacked 
objects used for the SED fitting are not exactly the 
same as those presented in M.l) but the stacking is per- 
formed only for objects with Spitzer/IRAC 3.6,4.5,5.8, 
and S.Ofim photometry from the Spzteer legacy survey of 
the UDS field (SpUDS; Spitzer Proposal ID 40021; PI: 
J. Dunlop) so that stellar population parameters be well 
constrainecFI. After removing objects with any confu- 
sion from neighboring objects in the IRAC images by 
eye, we stack 304 and 55 LAE candidates in the NB118 
sub-region and the NB209 sub-region, respectively. 

The procedure of t he SED fitting is the same as that 
of lOno et al.l ()2010b[ ). except for a fixed redshift (z = 
2.18). We use the stellar popu lation synthesis model 
GALAXEV (IBruzualfc Charlotl[2003') for stehar SEDs 
and include nebular emission (jSchaerer fc de Barros 
|2009| ). A Salpeter initial mass function flMF: [Salpetei 



We demand only IRAC coverage, not IRAC detection. 
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Fig. 9. — Results of SED fitting for tiie stacked objects in the NB118 sub-region (left) and the NB209 sub-region (right). The filled squares 
show the observed flux densities used for the fitting {B, V, R, i' , z' , J, H, K, [3.6], [4.5], [5.8], and [8.0]), while the open squares indicate those 
not used for the fitting {u* and NB387). The blue lines show the best-fit model spectra, and the blue crosses correspond to the best- fit 
flux densities. 



TABLE 6 

Broadband Photometry of the stacked LAEs for SED fitting 



sample B V R i' z' J H K [3.6] [4.5] [5.8] [8.0] 

NB118 sub-region 25.96 25.93 25.95 25.85 25.82 25.44 25.12 25.30 25.38 25.56 26.90 99.99 
(29.31) (29.28) (29.29) (29.20) (29.13) (28.60) (28.06) (28.41) (28.56) (28.08) (26.26) (26.00) 

NB209 sub-region 25.87 25.80 25.75 25.60 25.48 25.09 24.72 24.96 25.13 25.11 25.51 26.19 
(29.20) (29.11) (29.07) (28.93) (28.73) (28.06) (27.38) (27.86) (27.62) (27.28) (25.35) (25.04) 

Note. — Broadband photometry of the staeked objects whose stackings are performed for objects with IRAC coverage (see S[5j. All magnitudes 
are total magnitudes. 99.99 mag means negative flux densities. Magnitudes in parentheses arc Icr uncertainties adopted in SED fitting. 



TABLE 7 

Physical Properties of the stacked ob,iects from SED fitting 



sample Z E{B - V)i, Age Xr 
[Zq] [10'Mq] [nmg\ [Myr] 

NB118 sub-region O.rtg.l 0.2 2.88^0 13 0.2lt|j;[]4 ^'^-^t^f 1.678 

NB209 sub-region 0.8^0° 0.2 4.79lo;gi 0.27+'^,-°l 8.32^^;^^ 1.661 

Note. — Derived physical properties and their Itr errors of the staeked objects from SED fitting. Z is 
fixed to 0.2Zq. The degree of freedom is 8. The results may suffer from additional errors due to degeneracy 
with age or stellar metallicity, and due to possible systematic uncertainties. 



|1955() is assumed. We choose constant star formation 
history and the stellar metallicity Z = 0.2 Zq since pre- 
vious studies have shown that most LAEs are young, and 
that constant star formation history and subsolar stel- 
lar metallicities are reasonable assumptions. Since we 
include nebular emission in the fitting, varying star for- 
mation history a nd stellar metallicit y makes the fitting 
too complicated. IGuaita et al.l ()2011l ) performed SED fit- 
ting to 2 ~ 2.1 stacked LAEs with three scenarios of star 
formation histories, exponentially increasing, decreasing, 
and constant star formation, to find equally good fits to 
the data. They also note that among the free parame- 
ters SED fitting can relatively well constrain stellar mass 
and dust extinction, which are of particular interest in 
this paper. We also note that the assumption of subsolar 
metallicity is consistent with the gas-phase metallicity 



derived for our stacked LAE (see ij6.2l) . For dust ex- 
tincti on, we use Calzetti's extinction law (jCalzetti et al.l 
120001) on the assumptio i i oi E( B-V)^^ = E{B-V)^ as 
proposed bv lErb et al.l (|2006b[ ). IGIVI att enuatio n is cal- 
culated using the prescription given bv iMadaul ([1995). 
We do not use u* and NB387 data since they are sig- 
nificantly contaminated by strong Lya emission. Table 
[B] summarizes the broadband photometry of the stacked 
objects that are used for SED fitting. The uncertain- 
ties in optical and NIR band photometry listed in Table 
[B] contain photometric and two systematic errors. The 
photometric errors are estimated in the same manner as 
explained in ij4.1.1[ and the two systematic errors, asso- 
ciated with aperture correction and zero point, are given 
as follows; the errors in aperture correction are estimated 
to be 0.01 — 0.03 mag using point sources, while the zero- 
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point uncertainties are assumed t o be 0.05 mag for ever 
optical and NIR band (inferred bv lFurusawa et al] (|200< 
for optical bands). The photometric and two systematic 
errors are added in quadrature. For the IRAC channels, 
we include only photometric errors and do not include 
any systematic errors since they are unknown and the 
photometric errors are large enough to be dominant un- 
certainties. Inclusion of systematic errors to the IRAC 
channels will not change the results of the SED fitting 
significantly. More details of the SED fitting will be re- 
ported elsewhere (Y. Ono et al. in preparation). 
We obtain = 2.88lg;t| x 10^ M©, E{B - V) = 

0.2ll°j;^, and age = 1.261^;^^ x 10^ yr for the stacked 
object in the NB118 sub-region, with a reduced chi- 
squares of = 1-68, and = 4.79j:^:8i x 10^ Mq, 



that the [Nil] /Ha ratio varies with metallicity; indeed, 
this ratio is used as a metallici ty indicator of galax- 
ies, called the N2 index (e.g., iPettini fc Pagell 120041: 
iMaiolino et"alll2008t ) . In TableH we derive L(Ha+ [N ii] ) 
to be 8.01^ 2 X lO^'^ergs"^. We use the metallicity es- 
timat ed in g6j to infer £ ([Nii])/L(Ha) = 5.7til ^ 
10-2 (jMaiolino et al.ll2008( ). The Ha luminosity is then 
^V(l+5.7tl]xl0-2^ 



8.0;iJxlO 



from which we obtain SFR = M© yr 
With E{B - V) = 0.27: 



7.6l:i:2X 10^1 erg s-i. 



-1 



-0.01 
0.03 



combined with 



iCalzetti et al.l (|2000 r)'s extinction law, we obtain the 



E{B-V) = 0.27+'(j:^^, and age = 8.32+;;°° x 10^^ yr for the for typical z - 2 LAE s of L(Lva ) > 1 x lO'*^ erg s 



dust-corrected Ha luminosity to be 1.7_q;3 x 10 ergs , 
which is translated into SFR = 14^^ Mq yr"!. This SFR 
is the first unbiased SFR estimate from Ha luminosity 

Al- 



stacked object in the NB209 sub-region, with Xr = 1-66. 
Table [7] summarizes the physical properties derived from 
the SED fitting, and Figure IH] shows the best-fit model 
spectra with the observed flux densities. The errors in 
the best-fit parameters correspond to the la confidence 
interval (A^^ < 1) for each parameters. The relatively 
small errors of the results are due to the small uncer- 
tainties adopted in the SED fitting. We note that the 
results from the SED fitting may suffer from additional 
errors due to degeneracy with age or stellar metallicity, 
and due to other possible systematic uncertainties such 
as the position matching, and in the assumption of star 
formation history and stellar metallicity. The discussion 
will be expanded more in Y. Ono et al. (in preparation). 

Although the dust extinction of our LAEs is rather 
high compared with those of high er-z LAEs {E{B — 
V) = 0.00 - 0.07 at z ~ 3; e.g., iNilsson et ahl 120071 : 
lOno et al.|[2010aD . such an increasing trend of E{B — V) 
in LAEs down to z ^ 2 is also seen in other studies 
(INilsson et al.l 120111: [Guaita et al.ll201lD . INilsson et all 
pout ) performed SED fitting to 2 ~ 2.3 LAEs to find 
that Av {E{B - V)) varies over 0.0 - 2.5 (0.00 - 0.61) 
with an average of Ay = 0.6 {E{B - V) = 0.15) and 
Av = 1.5 {E{B -V) = 0. 37) for old aiid youn g popula- 
tion models, respectively. iGuaita et al.l ()201lD obtained 
E{B-V) = 0.22ta°t stacked z - 2.1 LAEs assuming 
constant star formation history. These findings combined 
with our result may indicate a strong evolution of dust 
extinction in LAEs from z > 3 to z ~ 2. 

6. RESULTS AND DISCUSSION 

6.1. Star Formation Rate 

6.1.1. Deriving SFR from the Ha Luminosity 

The Ha luminosity is believed to be the most reliable 
SFR indicator of galaxies among those based on the rest- 
frame UV and optical spectral features. Indeed it is pro- 
portional to the birth rate of massive stars as well as be- 
ing relatively insensitive to dust extinction as compared 
with UV-continuum. We measure the SFR of the stacked 
object in the NB2 09 sub-region fro m its Ha luminosity 
using the relation (|Kennicuttlll998l ): 



SFR[Mo yr"^] = 7.9 x ^-''^^(Ha) erg s 



(6) 



Before applying this relation, however, we have to 
subtract the contribution from [N ii]AA6584, 6548 lines 
from the observed Ha+ [N 11] luminosity. It is known 



though iFinkelstein et^TI ()2011|) have derived SFRs for 
two z ~ 2.3 — 2.5 LAEs from Ha spectroscopy, both ob- 
jects have L(Lya) ^ 2 x lO''^ ergs"^, which is about one 
order of magnitude brighter than typical L ya luminosi- 
ties o f LAEs from narrow-band surveys (e.g.. lOuchi et al.] 
[20i)8), du e to pre-select i on of bright objects for NIR spec- 
troscopy. iHaves et al.l ()2010l ) have measured L(Ha) for 
six LAEs with L(Lya) = (0.3 - 4.5) x lO*^ crgs^^ for 
which Ha emission is detected. In contrast to these stud- 
ies, our study is based on a large number {N = 105) of 
purely Lya-selected galaxies. 

SED fi'tting of the stacked object in the NB209 sub- 
region gives a stellar mass of 5 x 10^ M©. There- 
fore, the stacked object has the SFR and the stellar 
mass of the same order of magnitudes of t hose of LAEs 
studi ed earlier from SED fitting (z ~ 2: IGu aita et al.l 
20T1I z ~ 3: pawi s_cr ct al. 2006, 20p iLaiet al. 200^ 
Nilsson et ahl 120071: lOno et all [2010^ Some studies 



have derived SFRs of LAEs from (dust -uncorrected) UV- 
continuum (e.g. , GronwaU ct al. 200 f: lOuchi et al.ll2008L 



I2M0I : INilsson et al.ij200a . ,20 lit iGuaita et al.ll2010f ) which 
are roughly in the range 1 — lOMQyr"^, comparable to 
our dust-uncorrected SFR. 

We also estimate in a similar manner the SFRs 
of the (Ha-F[Nii])-detected objects, NB387-C-22326, 
NB387-W-04041, and NB387-W-04492 to be 23i^, 51^^, 

and 1313 ^0yr"\ respectiveljE3- These SFRs are 
much larger than that of a typical LAE at z: ^ 2, 
but similar to some bright LAEs at sim ilar redshifts 
(INilsson et al.ll2011l: IMnkelstein et al.ll201lD . This is be- 
cause our NIR narrowband images are relatively shal- 
low, and only bright, massive LA Es with high SFR 
are likely to be detected. Indeed, INilsson et aLl 's and 
IFinkelstein et ahl 's LAEs are much more massive than 
a typical LAE (> 10^° Mp) wit h some exceptions like 
HPS256 (IFinkelstein et al.l 1201 ID whose mass can be as 
smaU as 6 x 10® Mq. Although the (Ha-|- [N 11] )-detected 
objects lack Spitzer/IRAC data and their stellar masses 
are poorly constrained by SED fitting, rough estimate 

22 For NB387-C-22326, we assume L(Ha+[Nll]) = L(Ha), since 
its [0 11] is not detected and therefore its metallicity is low enough 
that the contribution of [Nil] to L(Ha+[Nll]) is considered to be 
neghgible. For NB387-W-04492, we adopt solar metalli city fr om 
the empirical [O ll]/(Ha-|-[N ll])-metallicity relation (see ^6.2. ip to 
correct for [Nil] emission. For NB387-W-04041, the metallicity is 
estimated to be 12 -f log(0/H) ^ 8.1 with Equation All SFRs 
are calculated assuming dust free. 
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This is consistent with the SFR de- 
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Fig. 10. — Specific SFR (sSFR) vs. stellar mass for various 
types of galaxies at 2 ~ 2. The red large filled circle indicates 
our result after correction for dust extinction of E{B — V) = 0.27, 
while the red horizontal bar indicates the sSFR on the assump- 
tion of E( B — V) = 0. The tria ngles are LAEs with Ha mea- 
surements IIFinke lstcin ct al. 20111). The squares arc UV-solccted 
galaxies (open: [Erb et al. 200614 filled: [Reddv ct al. 2006) and 
the pentagons arc sBzK galaxies (open: IHavashi et al.ii2009 . filled: 
lYoshikawa et al. 2010). All sSFRs are derived from Hq lumi- 
nosities after dust correction. A Salpeter IMF is assumed for 
all objects. The dashed lines correspond to constant SFRs of 
0.1, 1, 10, 100, 1OOOM0 yr-i. The gray filled circles with error- 
bars are z = 2 — 3 stacked LAEs whose SFRs are estimated from 
SED fitting JGuaita et al.ll2011|; INilsson e t"alll2007l : IGawiser eFall 
l200g. 12007: Lai ct al. 2008; Ono ct al. 2010"^), and the gray crosses 
are z = 2.3 LAEs w hose SFRs arc estimated from UV-continuum 
UNilsson et al.llMIIi '). 



is possible from their if-band magnitudes. iDaddi et al.l 
([2004), for instance, derived a relation between K mag- 
nitude and mass for z ^ 2 BzK galaxies. The K- 
band total magnitudes of NB387-C-22326, NB387-W- 
04041, and NB387-W-04492 are 23.28, 23.02, and 22.39, 
which correspond to a stellar mass of (2 — 3) x 10^°, 
(2 -4) X 10^°, an d (4-7) x 10^° Mq, respectively, follow- 
ing [DaddT^FaLl's simple relation. Although this is just a 
rough estimate and the relation itself involves a ce rtain 
uncertainty (cr(AlogM*) ~ 0.2: IDaddi et al.ll2Q04[ ). the 
individually (Ha-f [N ii])-detected objects are likely to be 
massive LAEs. 

6.1.2. Denying SFR from the [Oil] Lumtnosity 

The [O ii] luminosi ty is also known to b e a useful 
SFR indicator (e.g., [G allagher et al.' '1989"; 'Kennicutt' 
[1991 |1998; K ewlev et a l. 2004; Moustakas ct al. 2006) 
and frequently used to derive SFRs of galaxies at high- 
z, where Ha cannot be accessed from the ground (e.g . , 
iTeplitz et al'][2Q03l : lHopkinsll200l: iTakahashi etaI[[2Q071 ) . 
The SFR is derived from the [O ii] luminosity using the 
relation ([Kennicutt] 1 19981 ): 

SFR[Mq yr-1] = (1.4 ± 0.4) x IQ-'^^ L{[011]) erg s"\7) 

In Table [1 we derive i([Oii]) in the NB209 sub-region 



to be 5.3: 



\\ X IC*^ ergs ^, which corresponds to SFR 



rived from the Ha luminosity when dust free is as- 
sumed. However, if we take into account the dust ex- 
tinction of E{B -V) = 0.27t[] !]^, we find that the 
SFR derived from the [O ii] luminosity could be more 
than twice the SFR derived from the Ha luminosity, and 
they are not consistent within la errors: SFRcor(Ha) 
= Mt^Moyr-i, and SFReor([Oii]) = i2tl^MQyr-\ 
Such differenc es have also been re ported for present-day 
galaxies f e.g.. [Gilbank et al.|[2010D and for high- z galax- 
ies (e.g., [ChariotetaDlloOllSei^ The 
difference can be due to the complicated origin of the 
[O ii] luminosity. The [O ii] luminosity is not directory 
proportional to the ionizing luminosity, and also depends 
on the chemical a bundance and excitation state of the 
ionized gas (e.g [Kennicutti [1998[: [Kewlev et all [2004[: 
[Moustakas et alj [2OO60. Indeed, the ratio of [Oii]/H/? 
can b e used as a metallicity indicator (e.g.. [Nagao et alJ 
[20061 : see also ^6.2.ip . Alternatively, the uncertainty of 
the dust correction may cause the difference of the SFRs. 
The relatively large dust correction for the [O ii] luminos- 
ity makes it difficu lt to derive the SFR accurately (e.g., 
[Jansen et all[200l . The difference of the SFRs of our 
result may be also due to the overestimates of the dust 
extinction inferred from the SED fitting. In any case, we 
use the SFR derived from the Ha luminosity, which is 
more directory proportional to the SFR and less affected 
by the dust attenuation or other factors (e.g., metallic- 
ity), in the following analysis. 

6.1.3. Mass-sSFR Relation 

We plot the stacked object on the specific SFR (sSFR 
= SFR/Af^) vs. stellar mass plane in Figure (TUI Al- 
though sBzK and UV-selected galaxies obey a simple 
scaling relation between the sSFR and the stellar mass, 
our stacked LAE is located below an extrapolation of this 
relation toward lower-mass, even though the dust extinc- 
tion of E{B — V) = 0.27 is taken into account in the SFR. 
We note that the other LAEs whose SFR and stellar mass 
are derived by SED fitting (gray symbols) are also dis- 
tributed below the extrapolation. If LAEs are typical 
of low-mass galaxies, this figure indicates that low-mass 
galaxies with M < 10^^^" A/0 have lower star formation 
efficiencies than extrapolated from more massive galax- 
ies. This trend, if real, is qualitatively consistent with 
galaxy formation models which predict that less massive 
systems have low star formation efficiencies owing to var- 
ious mechanisms such as feedback from supernovae (see, 
e.g.. lBenson.,20101 for a review). 

6.2. Metallicity 

6.2.1. Constraining the Metallicity 

An accurate estimate of gas phase metallicity re- 
quires knowledge of the electron temperature which is 
provided by comparing auroral lines to nebular emis- 
sion lines (e.g., [O iii]A4363/A5007; direct method; 
[Lee et ahl [2004) . However, auroral lines are generally 
weak and it is difficult to observe them in distant galax- 
ies. A number of empirical relations between the ra- 
tio of nebular lines and metallicity have been proposed 
to measure gas-phase metallicities of distant galaxies. 
Among them, we combine the [Nli]A6584/Ha indicator 
{N2 index: e.g., [Pettini fc Pagell \2004 [Maiolino et"all 
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Fig. 11.— (bottom) Relation between [O ll] A3727/(Ha+ 
[N ll] AA6584, 6548) and gas-phase metallicity. The red solid line 
corresponds to the empirical relation given by Equation ^ com- 
bined with Equation JgJ. The gray diamonds represent low- 
metallicity galaxies compiled by Nagao et al. ( 20j3) whose metal- 
licities a re inferr ed through the direct Te method and then used 
in Maiolmo et ah, C200S ) to fit the polynomials in Equation (|9]l 
in the range 12 + logZ < 8.3; they are distributed around the 
best-fit line with a large scatter (A l og Z ~ 0.2). The gray dots 
are SDSS galaxies IIKewlev fc Dopit"a|[2002i) used for the fit in the 
range of 12 + logZ > 8.3 whose metallicities are inferred through 
photoionization models. The dashed and long-dashed lines cor- 
respond to the relations derived from a combination of stellar 
population synthesis and photoionization models with an ioniza- 
tion p arameter of q = 1.5 X 10^ and 1.0 X lO^cms"^, respec- 
tively I IKewlev &c Dopit"al 120021) . We ignore the range Z > Zq 
in this paper, beca use L AEs are found to be metal-poor in pre- 
vious studies (see il6.2.1| l. The blue shaded area is the range of 
[Oll]A3727/(Ha-|-[Nll]AA6584, 6548) for the stacked object assum- 
ing E{B — V) = 0. (top) Residuals from the best-fit polynomials 
for the low-metallicity galaxies (left) and their histogram (right). 
The two dotted lines show the rms of the residuals for A log Z > 
objects (rms= 0.22) and AlogZ < objects (0.26), respectively. 



20081) and the [Oii1/H/3 indicator (e.g., INagao et al] 
200a iMaiolino et"al[|2008[ ) to derive the metalhcity of 



our object. Even for the same metallicity index, how- 
ever, calibrations are often different among the authors, 
leading to di fferent metallicity measures as discussed in 
INagao et all 12006). As for the N2 inde x INagao et al.l 
(l2006D's calibration and IPettini fc Pagel (I200j)'s agree 
well (AZ < 0.2 dex) with each other over 7.7 < 12 -|- 
log(0/H) < 8.5, while out of this range a large differ- 
ence is seen probably due to the lack of objects for cal- 
ibration in l^ttini & Pagcl (2004) sample. We use the 
calibrat ion bv j Mai olino ct ah, ((2008t ). which is an up- 
date of INagao et al . (2006), based on a \aw-Z sample 
(7.7 < 12-Hlog(0/H) < 8.3). 

We cannot, however, apply these indicators di- 
rectly to our object, because we have only [O ii] and 
Ha-|-[Nii] fluxes. Instead, we make use the fact that 
[Oll]/(Ha-|-[Nll]) ratio can be expressed as a combina- 
tion of the two indicators: 

\Oll] 



2.85 (1 -I- 1.33 X [NII]A6584/Ha) ' 



(8) 



where we assume that all lines are dust free and 
adopt an intrinsic Ha/H/3 ratio of 2.85 and an in trinsic 
[Nii]A6584/A6548 ratio of 3.0 (|Osterbrockl[T989h . The 
effect of dust extinction will be discussed later. The 
metallicity dependence of the two indicators is empiri- 
cally approximated by the polynomial: 



log R~Cq-{-CiX + C2x'^ 



C3X 



(9) 



where R = [O ii]/H/3 or [N 11] /Ha, x is the metallicity rel- 
ative to solar (i.e., x = log(Z/Z0) = 12+log(0/H)-8.69, 
lAllende Prieto et al.l 1200 li) . and the coefficient s cn t o 
C4 are taken from Table 4 of IMaiolino et al.l (|2008| ). 
Thus, [O ii]/(Hq;-|-[Nii]) ratio is expressed as a func- 
tion of metallicity, as shown in Figure [TT] by the solid 
curve. In this figure, the gray diamonds and gray dots 
indicate, respectivel y, local low-nietalli city galaxies with 
12+log(0/H) < 8.3 (|Nagao et alMM\ and SPSS galax- 
ies with 12 + los(0/H ) > 8.3 (jTremonti eraLll200l 
iKewley fc Dopital 120021 ) used to calibrate Equation dH); 
the metallicitie s of the former are measured with the di- 
rect Te method (jNagao et al.ll2006l ) , while those of the lat- 
ter are derived by applying photoionization models to the 
most prominent optical emission lines ([On], II/3, [Oiii], 
Ha, [Nil], [Sii]). The ratio has a peak at around solar 
metallicity and for any given value of [O ii]/(Hq;-|-[N 11]) 
(except for the peak), there are two solutions of metal- 
licity, one being subsolar and the other being super- 
solar. The blue shaded region in Figure [TT] corre- 
sponds to the observed [O ii]/(Ha+[N ii]AA6584, 6548) 
ratio including the la photometric error: 
Two metallicity ranges are found to meet the observa- 



8.21 



+0.10 



and 8.94 



+0.05 



The 



(Ha+ [NII]AA6584,6548) 



tion: 12 + log(0/H) - o..._o.ii ^..^ "---o.os; 
latter range, however, appears to be unlikely, since it 
is not consistent with recent spectroscopic observations 
that L AEs have much lower metallicities than t he solar 
value (jFinkelstein et al.l 1201 It iCowie et al.ll2011h . Mod- 
est dust extinction inferred from SEP fit t ing a l so fa- 
vors low m etallicit ies (e.g., 'Gawisc r et all 120 06. 2007t 
iPirzkal et al. 2007: iLai et al.i 2008: Quo et ahfe oiOa). If 
we rule out the supersolar solution, the metallicity of our 
object is estimated to be 12 + log(0/H) = 8.2ll°;J°, or 

z/Zq = o.ssIq'q^. 

This metallicity estimate, however, involves two sys- 
tematic uncertainties. First is the uncertainty in the 
polynomial fit given by Equation ([9]) for low-metallicity 
galaxies. The gray diamonds plotted in Figure [TT] 
represent galaxies w ith 12 -I- log(0/H) < 8.3 used in 
IMaiolino et al.l ([2008 1 to fit the polynomials. The galax- 
ies are distributed around the solid line (as expected), 
but with a large scatter of AlogZ = -^0.22/- 0.26 (Ict). 
This may suggest that there are large errors in mea- 
surements of line ratios and/or metallicities, since these 
measurements collected from the literature are b ased on 
differ ent methods to measure line ratios (Nagao et al.l 
|2006| ). Alternatively, such a large scatter may be in- 
trinsic. Pashed lines in Figure [TT1 represent relations of 
metallicity and [O ii]/(Ha-|-[N 11]) from a combination of 
stellar population synthesis and photo-ionization models 
with a set of ionization parameter {q; Kewley & PopiTal 
l2002f ). From Figure [TT] the large scatter could be ex- 
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TABLE 8 

Coefficients for metallicity indicators in Equation ([9j 



Flux ratio (log R) cq c\ C2 c-j C4 

log([Nii]A6584/Ha)(i) -0.7732 1.2357 -0.2811 -0.7201 -0.3330 

log([Oii]A3727/H/3)(i) 0.5603 0.0450 -1.8017 -1.8434 -0.6549 

log([Oii]A3727/(Ho+[Nii]AA6584, 6548))(2) 0.0164 -0.1673 -1.9484 -1.8158 -0.6148 

^ ^ The coefficients are taken from Table 4 of lMaiolino et all 1200SD . 

(2) The coefficients are determined by a best fit polynomial of the combined equation described in Equation ( fSl . 




from R23 from N2 



12 + log(0///) 

Fig. 12. — Comparison of the metallicity calculated from our 
[O ll]/(Hci+[N 11]) indicator with those from two frequently used in- 
dicators, the R23 index (left panel) and the N2 index (right panel), 
for z ~ 2 Icnsed galaxies taken from the literature (Rig bv et al. 
2011 with triangle; Richard et al. 2011 with squares; Hainli ne et al. 
2009 with pentagons). At each panel, the dashed line is the 
line of equality, and the dotted lines show the la errors associ- 
ated with our indic ator in cluding the calibration errors (A log Z = 
-1-0.22/ — 0.2 6; see i|6.2.1l l. Our indicator has in principle two so- 
lutions (see i|6.2.1l l. and for galaxies with two solutions over their 
1(T errors both metallicities are plotted; solutions with the smaller 
differences from Zfl23 or Zjv2 in black, while those with the larger 
differences in white. 

plained by a diversity of ionization parameters {q 
(1 — 10) X 10^ cm s"-'^ ). Note that the em pirical line can be 
reproduced by the iKewlev fc Dopital 's photo-ionization 
model with q ~ (4 — 8) x 10''cms~^. In any case, 
we estimate that the calibration error in our metallic- 
ity estimate due to the uncertainty in the polynomial 
fit is AlogZ = +0.22/ - 0.26 (Ict). In Figure [H we 
compare metallicities derived from our original indicator 
with those from more commonly used indicators: the 
R23 index ff[Oii1A3 727-HfOiii]AA4959,5007)/H^: e.g., 
iTremonti et al.ll2004l and references therein) and the N2 
index to check robustness of the indicator. We com- 
pare the indicator s using spectroscopic data of z ^ 
2 lensed galaxies (iHainline et al.l 120091 : iRichard et al.l 
[20TI iRigbv et all 1201 ID . From the comparison, we see 
that metallicities from our indicator are roughly consis- 
tent with those from other indicators within A log Z — 
-fO.22/- 0.26 (Ifj). 

Second is the uncertainty in dust extinction. While 
we have assumed E{B — V^) = to derive Equation (|8]), 
the SED fit suggests that the stacked object may have 
E{B - V) up to ~ 0.27 (see We find that adopting 
-V)= 0.27 instead of E{B - V) ^ increases the 
metallicity estimate by AlogZ ~ 0.3. 

When these two systematic errors are taken into ac- 
count, the metallicity range of our object is estimated to 
be 12 + log(0/H) = 8.2lt°;J° (random) (calib) = 
8.2ltolj, or Z/Zq = O.SSto fe^ assuming E{B^V) = 0. 



If E{B — V) = 0.27 is adopted, the metallicity increases 
to 12 + log(0/H) ~ 8.5toit, or Z/Zq = O.Ttoi- Re- 
call, however, that our method logically permits the pos- 
sibility of a supersolar metallicity and that in order to 
rule out such a possibility, we will require independent 
data favoring metal-poor LAEs. In this sense, the upper 
limit of the metallicity obtained above is not as strict as 
the lower limit. Considering this, we adopt a conserva- 
tive conclusion that the metallicity of the stacked object 
is no less than 0.17 {la), or 0.09 {2a). This 2a 
lower limit is obtained in the following manner. The ob- 
served [O ii]/(Ha-f [N 11]) ratio and its 2a negative pho- 
tometric error are 0.66 and —0.23, respectively, thus the 
2a lower limit of the observed [O ii]/(IlQ;-t-[N 11]) ratio is 
0.43, which corresponds to 12-f log(0/H) = 7.96. There- 
fore, the 2a negative statistical error in terms of metal- 
licity is (AlogZ)caiib = -(8.21 - 7.96) = -0.25. On 
the other hand, the 2a negative systematic error is just 
twice the la value (—0.26), hence (AlogZ)sys = —0.52 
{2a). Therefore, we obtain the 2a negative error in 
terms of metallicity as (AlogZ) = -^0.252 -t- 0.522 ^ 
—0.58, and the 2a lower limit of the metallicity as 
12 + log(0/H) = 8.21 - 058 = 7.63, corresponding 
to Z = 0.09^0. We calculate the 3a, 4tT, and 5cr 
lower limits of metallicity in the same manner to obtain 
12 + log(0/H) = 7.32, 6.97, and 6.58, or Z = 0.04 Zq, 
0.02 Zq, and 0.008 Z©, respectively. Although a recent 
spectroscopic study has placed a weak upper limit on 
the metallicity of brig ht LAEs (Z < 0.56 Z© (2cr0; 
iFinkelstein et al.]l2011h . this is the firs t lower limit on the 
metallicity for typical LAEs at high-z. iScannapieco et aTl 
(|2003[) proposed that LAEs may b e extre mely metal poor 
primordial galaxies, and iSchaereii (|2003D also found that 
large Lya EWs can be observed from extremely metal 
poor galaxies (Z/Z© < 10~^). However, our firm lower 
hmits, e.g., Z > 2 x 10~'^ Zq at the 4tT level, do not 
support this idea at least for z = 2.2 LAEs. For higher 
redshift (> 3) LAEs, on the other hand, their metallici- 
ties may be lower than what we derive for z = 2.2 LAEs, 
because, as seen in SJSI the amount of dust extinction for 
LAEs decreases with redshift from 2 '-^ 2 to z > 3. 

6.2.2. Mass- Metallicity Relation 

Figure [13] shows the observed mass-metallicity {M — Z) 
relation of star-form ing galaxies at loc al and z ^ 2.2 com - 
piled bv lMaiolino et al. (2008): Kcwl ev fc ElhsonI (pOOl 
for z ^ 0.07. and [Erb et al.. (,2006ail for z ~ 2.2. For 
z ^ 2.2, the relation is derived based on stacked galax- 
ies with M > 3 X 10® M©. The line below the observed 

W e have recalibrated the metallicity given in Finkelstei n et aLl 
II2OIII ) using Equation The original metallicity upper limit is 
Z < 0.41 Zq (2o-). 
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Fig. 13.— Mass-Metallicity (M - Z) relation of UV-selected 
galaxies and LAEs at z ~ 2. The red circle with an arrow 
represents our stacked LAE. The red ticks show the lower lim- 
its of metallicity at 1 and 2ct levels as labeled next to the ticks. 
The open triangles show two LAEs at 2 ~ 2.3 and 2.5 from 



IFin kelstein et aTl I I2011I) whose upper limits denote the 2a confi- 
dence limits. The curves indicate the M — Z rel ation observed 
at 2 ~ 0.07 (long-dashed; 'Ke wlev fc EiiisonI 120081 ') and z ~ 2.2 
(dashed: | Erb ct al. 2006a) (best-fit functions are determined by 
IMaiolino eFal., 1200S') ). The filled squares are stacked data point; 
of 2 ~ 2.2 UV-selected galaxies in six bins in stellar mass(Erb ct 
| 2006al ). and the open square indicates BX418 {z = 2.3; ,Erb et 



1201 0). which is considered to be the most metal-poor UV-selected 
galaxy. The gray diamond area shows the median stellar mass 
and metallicity for z = 0.195 — 0.44 LAEs and their Itr ranges 
llUowie e t al. 20lJ). Al l data have been reca librated to have the 
same metallicity scale UMaiolino et al.l 120081 ) and IMF IjSalpeterl 
119551 ) so that all results can be directly compared. 

limit of stellar mass is therefore a smooth extrapolation 
of the best-fit function. The M — Z relation is found to 
evolve with redshift in the sense that the metallicity at a 
given stell ar mass decreases with increasing redshift u p 
to z ~ 3.5 (jMaiolino et al.ll2008t iMannucci eralll2009f ). 

The red circle with an upward arrow represents our 
stacked object in the NB209 sub-region which has a stel- 
lar mass of 5 X 10^ M© and a metalhcity of Z > 0.09 Z© 
(2it) . The large circle represents the central value, which 
is derived assuming E{B — V) = and the two ticks 
stand for the la and 2a lower limits. Note that our 
study is the first to place a lower limit to the metallic- 
ity of z ^ 2 galaxies with stellar masses below lO^M©, 
well below the lower-mass limit of the previous studies. 
The two triangle s ymbols correspond to the two LAEs at 
z - 2.3 and 2.5 bv IFinkelstein et al.l (|2011| ): for both ob- 
jects the metallicity estimate is an upper limit. The less 
massive one o f the two has a st ellar mass comparable to 
our object. If ' Finkelstein et all observed a typical LAE, 
then the combination of their results with ours suggests 
that LAEs with < lO^M© at z ~ 2 have metallicities in 
the range 0.09 < Z/Zq < 0.56 with the 95% confidence 
level. 

Comparison of our object with the M — Z relation 
at z = 2.2 reveals that our object is consistent with a 
smooth extrapolation toward lower masses of the M — Z 
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Fig. 14. — Fund amental M ~ Z relation proposed by 
IMannucci et al.1 I I2010I ) . The dashe d curve indicates the best-fit 
relation for 2 ~ SDSS galaxies l IMannucci et al.|[20Tol . [2011). 
whose typical distribution ranges are shown as the gray shaded 
(for those with > lO^-^ Mg, 

chabrier) the blue shaded 

{Mi, < 10^'^ Mq chabrior) regions. The black circles are z ~ 2.2 
UV-selected galaxies compiled bv lMannucci et al.J l|201 0|), the blac k 
squares are 1.5 < z < 2.5 lensed galaxies micWd etaIll20Tll ), 
and the ope n triangles show two LAEs at 2 ~ 2.3 and 2.5 from 



IFinkelstein et aTl I I2011I ). The red circle with an arrow is our re- 
sult; it is located near a smooth extrapolation of the fundamental 
M — Z relation toward lower stellar masses. 



relation at the 1.5a level. However, we also note that the 
central value of our object is 0.4 dex larger than the ex- 
trapolation. This offset might be real, since our metallic- 
ity estimate is a conservative lower limit. If real, there are 
two possible explanations of this offset. One is that the 
slope of the Al — Z relation may become shallower below 
the stellar-mass limit of UV-selected galaxies and our ob- 
ject is in fact on the relation, implying that UV-selected 
galaxies and LAE s obey a common M — Z relation. Al- 
though [ErbIeF|n] ()2010f ) have recently found an unred- 
dened, low-metallicity, and low-mass UV-selected galaxy 
(BX418) dX z — 2.3 (black open square) which seems 
consistent with th e conventional M — Z relation at 2; ~ 2 
(|Erb et al.ll2006a[) . it is only one object and its upper 
limit for UV-selected galaxies is still larger than our lower 
limit. The other possibility is that UV-selected galaxies 
and LAEs obey different M — Z relations and that the ex- 
trapolation is valid only for UV-selected galaxies. In this 
case, it is likely that LAEs have relatively high metal- 
licities for their masses, unlike the relation estimated for 
UV-selected galaxies. Although this appears to be in- 
consistent with the conventional picture that LAEs are 
the most metal-poor population, there are, in fact, few 
studies which compare metallicities of LAEs and other 
galaxies statistically at the same stellar mass. Thus, we 
cannot im mediately rule out this p ossibility. 

Recently, IMannucci et al.l (|2010f ) have found that the 
observed dispersion in the M — Z relation is correlated 
with the SFR in the sense that galaxies with lower SFRs 
have higher metallicities, and that star-forming galax- 
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ies at all redshifts below z ^ 2.2 obey a common, sin- 
gle M — Z— SFR relation. They referred to the re- 
lation as the fundamental M — Z relation. A simi- 
lar relation bet ween the three quantities is reported by 
ILala-Lopez et al . (2010,). a nd this trend is also realized 
for high-z lensed galaxies (jRichard et al.|[20TT[ ). Indeed, 
the SFR of our object is lower than that of UV-selected 
galaxies with similar mass, and most interestingly, we 
find that our object is located near a smooth extrapo- 
lation of the fundamental M — Z relation toward lower 
stellar masses (Figure [Til) . In Figure[T31 the dashed curve 
shows th e fundamental M — Z relation defined by SDSS 
galaxies (jMannucci et al.|[2010ll2011^ ■ and the gray and 
blue shaded areas show the typical distribution ranges 
of galaxies with Mi, > 10^-^ Afo^chabrier (|Mannuc ci et al.l 
[20lb) and < lO^ ^ Mo,chabrior (jMannucci et al.ll201lD . 
respectively. Therefore, the offset seen in the M — Z re- 
lation may be due to the relatively low SFR of a typical 
LAE, and the offset seen in the mass-sSFR relation may 
be due to the relatively high mctallicity for its mass com- 
pared to the value estimated from the M — Z relation of 
UV-selected galaxies. 

6.3. Lya Escape Fraction 

We infer the escape fraction of Lya photons (/^c") 
from the Lya and Ha luminosities. The Lya escape frac- 
tion is an important quantity of LAEs because it can be 
used to probe distribution of ISM in LAEs. Since Lya 
photons are resonantly scattered by neutral hydrogen 
(HI) gas in the ISM, f^/^" is strongly dependent on kine- 
matics and distribution of the ISM as well as the mctal- 
licity of the ISM. For example, f} ^-" wi l l be larger if some 
ISM i s outflowing (e.g. , Kunth et al.l 119981: lAtek et al.l 
I2008t iDiikstra fc Wvithe. ,2010, ). Clump y distributions 
of the ISM also makes & larger (e.g.. [Neufeldl [l99ll : 
[Hansen fc Oh| [2006': "F mkelstein et all[2008[r 

The Lya escape fraction is calculated by dividing the 
observed Lya luminosity by the intrinsic Lya luminos- 
ity produced in galaxy due to star formation. To obtain 
iint(Lya), most studies have used SFRs derived from 
SED fitting or UV continuum emission. However, in- 
stantaneous SFRs from SED fitting are model-dependent 
and can have erro rs as large as two orders of magnitude 
(|Ono et al.ll2010af ). and deriving SFRs from UV continua 
requires large correction for dust extinction. As men- 
tioned before, the SFRs are estimated most reliably from 
Ha luminos ities. 

Recently, IHaves et al.l (|2010f ) have compared luminos- 
ity functions of LAEs and Ha emitters (HAEs) at z — 2.2 
and found the volumetrically averaged escape fraction for 
star- forming galaxies to be ~ 5 %. This estimate is based 
on intrinsically different populations (LAEs and HAEs) 
and thus it is not clear whether or not LAEs typically 
have such low escape fractions. 

We estimate /^^" of LAEs at z ~ 2.2 by directly com- 
paring the Lya luminosity of the stacked object with its 
Ha luminosity. We calculate /^^" as: 



1042 



10« 



Lya ^ i^obs(Lya) ^ £obs(Lya) 
~ Ant(Lya) 8.7iint(Ha)' 



(10) 



where subscripts 'int' and 'obs' refer to the intrinsic and 
observed quant ities, respectively, a nd we assume Case B 
recombination ()Brocklehurstl[l97ll ). 
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Fig. 15. — (top) Lya escape fraction {fasc"') of LAEs against the 
observed Ly« luminosity. The red filled circle indicates our stacked 
LAE after correction for dust extinction of E{B — V) = 0.27, and 
the red horizontal bar is the case assuming E{B ~ V) = 0. The 
cyan circles are three individually Ho detected LAEs; for them we 
assume E{B — V) = 0, implying that their /<!^c" values are upper 
limits; if they have E{B — V) = 0.27 as in the case of the stacked 

object, decreases by 56% as indicated by the cyan arrow. 

The filled squares show ^ = 2.2 LAEs which arc also selected as 
Hcf emitters by a double narrowband survey (Hayes et al. 2010). 
The open triangles denote LAEs at z ~ 2.3 and 2.5 whose Lya and 
Ha luminosities are spectroscopically measured ( Finkelste in et al.l 
120111) . The dashed line shows fe^°'= 1. (bottom) Intrinsic Ha 
luminosity against the observed Lya luminosity. The meanings of 
the symbols and the line are the same as in the top panel. 



As listed in Table ffl Lobs (Lya) is 1.80t;^:EJ^ x lO*^ 
ergs-i. Lint(Ha) is derived to be l-lStoH x lO^^ergs"! 
by correcting Lobs(Ha) (7.55l};![5 x 10'*^ ergs~\ see Ta- 
bleg]) for dust extinction oiE{B-V) ^ 0.2?!°;°^. From 
Equation we find f^/^"^ 12^1 % for the stacked ob- 
ject. When dust free is assumed, /^^" can be as high 
as 27^4 %, which is about six times higher than that de- 
rived volumetrical ly for star- forming galaxies at z = 2.2 
(jHayes et al.ll2010l ). These values suggest that for a typ- 
ical LAE at z '--^ 2 a relatively large fraction of Lya pho- 
tons can escape the galaxies. Our value is also much 
higher than thos e estimated f or LBGs at z ~ 3 (e.g., 
~ 5% (median): iKorneil |201(][ ). but similar to that for 
LAEs at z = 2 - 3 (29% (median): iBlanc et all 1201 IL 
> 32%: IHaves et al.l [2010l using LAEs alone, > 14 %: 
IZheng et aLll2011l ). 

It is intuitively unreasonable that galaxies with E{B — 
V) = 0.27 c an escape as high as 12% of Lya pho- 
tons. Indeed, IKorneil (|2010D find that z ~ 3 LBGs with 
E{B-V) = 0.2-0.3 have f^^" of just a few percents. As 
described in ^ however, there is a possibility that dust 
does not always attenuate Lya emission, but rather plays 
an important role in escaping Lya photons from galaxies. 
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In order to quantify the effect of dust on flZ" , intro - 
duce the parameter q following iFinkelstein et ahl ()2008[ ). 
which is defined as q — T(Lya)/Ti2i6, where T(Lya) and 
T1216 are defined as e~"^^^^°'^ = iobs (Lya)/-£'int (Lya) and 
g-Ti2i6 = iQ-0Aki2i6E{B-v) ^j^j^ f,^^^^ being the extinc- 
tion coefficient at A = 1216 A. Small values {q < 1) mean 
that Lya photons suffer less attenuation by dust than 
UV-continuum photons, as expected for a clumpy distri- 
bution of the ISM (e.g., Ncufcld 1991) or spe cial ki nemat- 
ics of the ISM (e.g., outfiows; ,Kunth et aIlll998D . while 
large values {q > 1) means that Lya photons are more 
heavily attenuated by dust. 

The q parameter is expressed in terms of /^^" and 
E{B ~ V) as: 



LAEs. 



7. CONCLUSIONS 



iog(/e^r) 



0Aki2wEiB -V)' 



(11) 



Using fci2i6 = 11.98 (|Calzetti et al.l[200a we obtain q 



0.7 



+0.1 

-0.1- 



Similar results have been obtained from other 



studies. For example, iHaves et al.l (120101) obtain ed q ~ 
1 - 1.5 for z = 2.2 LAEs. and lBlanc et all (|2011[ ) found 
a median oi q ^ 0.99 for their z = 2 — 4 LAE sample. 

We note again that the E{B — V) value of our stacked 
object estimated by SED fitting may suffer additional 
errors due to possible systematic uncertainties (see 
so may the g-value. However, even when we take an 
extreme case of E{B — V) = 0.1, q increases only up to 
'--^ 1.5, implying that very large q (s> 1) are unlikely. Our 
result thus favors models of LAEs' ISM being outfiowing, 
or possessing a multi-phase clumpy distribution, or both, 
rather than in a homogeneous, static distribution. 

We also apply Equation pD)) to the individually de- 
tected objects, NB387-C-22326, NB387-W-04041, and 
NB387-W-04492 to obtain /Ly"= 6+1, 5^1, and 16t^ %, 
respectively. These values are derived from dust uncor- 
rected Ha luminosities, and hence are upper limits. Fig- 
ure [15] shows /^^" against the observed Lya luminos- 
ity for our LAEs (the red circle for the stacked object 
and the cyan circles for the individually detected objects) 
together with those LAEs taken from the literature for 
which Ha emission is individually detected: six LAEs at 
z ~ 2.2 whose H a luminosities are estimated from nar- 
rowband imaging (jHaves et al.ll2010D and two LAEs wit h 
Ha spectra at 2: ~ 2.3 and 2.5 ( Finkelstein et al.ll2011[ ). 
It is found that most of the LAEs with individual Ha 
detection have fcZ"^ 10%, with objects brighter in Lya 
luminosity tending to have higher escape fractions. It 
is interesting that our stacked object is roughly on this 
trend, since LAEs with individual Ha detection should 
be biased toward higher Ha luminosities. It may imply 
that most LAEs (irrespective of Ha detection) have rela- 
tively high /^^", at least higher than the vqlumet rically 
averaged star- forming galaxies (jHaves et al.ll2010l) . The 
trend of increasing f^J^" with Lya luminosity seems to be 
natural because galaxies with higher /^^" have brighter 
Lya luminosities and thus are more easily selected as 
LAEs. Finally, we note that recent studies have shown 
that a fraction of Lya is emitted from diffuse outer halos 
of galaxies (e.g., iSteidel et al.ll201ll ). We estimate the 
strength of Lya from the color u*— NB387 with a cer- 
tain aperture, therefore the escape fraction we derive is 
only the fraction escaping from the central region of the 



We have presented the results of the first detection 
of [0 11] and Ha emission from a typical Lya emitter 
(LAE) at z = 2.2 using a stacking analysis of a sample 
constructed from our Subaru/Suprime-Cam narrowband 
(NB387) survey in the Suharn/ XMM-Newton Deep Sur- 
vey field. The redshift z — 2.2 is unique, because [O 11] 
and Ha lines fall into NIR wavelengths where OH-airglow 
is very weak. We found 919 LAE candidates in this field. 
Follow-up spectroscopy was made for 30 candidates se- 
lected to cover wide ranges of NB387 magnitudes and 
u*— NB387 colors. Among the 13 out of 30 candidates 
with NB387 < 25, 10 were confirmed as z 2.2 LAEs. 
No emission line was detected for the remaining 20 with 
NB387 > 25 due to the lack of sensitivity. 

The near infrared observation was made to detect [O 11] 
and Ha lines by the NewHa Survey (Lee et al. in prepa- 
ration) with KPNO/NEWFIRM, using the narrowband 
fihcrs NB118 and NB209, respectively; 561 LAEs are lo- 
cated in the area covered by NB118 (NB118 sub-sample), 
among which 105 have also NB209 imaging (NB209 sub- 
sample). 

Only seventeen and seven candidates are individually 
detected in NB118 and NB209, respectively. However, a 
stacking analysis of a large number of undetected LAEs 
yielded statistically significant detection in both narrow- 
bands. We estimated the [On] and Ha+[Nii] EWs and 
fluxes of the stacked object using Monte Carlo simula- 
tions, and used the estimates to derive the SFR, gas 
phase metallicity, and Lya escape fraction (/^^") of 
a typical LAE. Our main results probed by the triple 
narrow-band survey are summarized as follows. 

• The Ha luminosity of the stacked object, after cor- 
rection for a contribution from [Nil] lines (esti- 
mated to be modest) to the Ha+[Nii] photometry 

+0.01 

03 



and for a dust extinction of E{B - V) 0.27Za, 



(derived from SED fitting), is l.Ttoi x 10''^ ergs-\ 
which corresponds to a SFR of 14^3 Mq yr~^. This 
is the first estimate of the SFR of a typical LAE at 
high-z based on the Ha luminosity. Adopting an 
stellar mass of 5 x 10^ A/0 derived from SED fit- 
ting, we plot the stacked object on the specific SFR 
(sSFR) vs. Mi, plane, and find that our stacked ob- 
ject is located below a simple extrapolation toward 
lower-masses of the observed sSFR - Mi, relation of 
z ^ 2 BzK and UV-selected galaxies. This trend 
is also evident for LAEs at similar rcdshifts whose 
SFRs are inferred by SED fitting. This indicates 
that low-mass galaxies with M < 1Q^~^^ Mq have 
lower star formation efficiencies than expected from 
massive galaxies. 

We use the hne ratio [0 11]/ (Ha+[N 11]) as a metal- 
licity indicator, and find that the metallicity of the 
stacked object is no less than 0.09 Zq at the 2a 
level. This is the first constraint on the metallic- 
ity of a typical LAE at high-z, and this relatively 
high lower- limit does not support, at least at z '-^ 2, 
the hypothesis that LAEs are extremely metal poor 
(Z < 2 X 10"^ Zq) galaxies at the 4cr level. We plot 
the stacked object on the mass-metallicity {M — Z) 
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plane, and find that the stacked object is not con- 
sistent with a simple extrapolation toward lower 
masses of the observed M — Z relation of z ~ 2 
UV-selected galaxies. Instead, our result seems to 
be consistent with the recently p roposed fundamen- 
tal M — Z relation ( Mannucci et al. 2010) for which 
the relatively low SFR of the stacked object is taken 
into account. 

• From the Lya and Ha luminosities, we found that 
the fig"" of the stacked object is 12^^ %, and can be 

as high as 27^4 %, much larger than those inferred 
for yolumetrically averaged star forming galaxies 
at z = 2.2 and Lyman-break galaxies at higher-z, 
but comparable to those of LAEs at z = 2 — 3. We 
compiled the f^/c"' data of LAEs with Ha emission 
from the literature, and found that most LAEs have 
relatively high flZ"iZ 10%), and that there are a 
mild trend that brighter LAEs have higher /^^". 
We also found a low q value for our object, q = 
0.7;to i- AH these findings indicate that LAEs have 
some unique mechanisms to efficiently emit Lya 
photons. 



The NB387 data used in this work are collected at the 
Subaru Telescope, which is operated by the National As- 
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